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USERS' 


MANUAL FOR COMPUTER PROGRAM FOR THREE-DIMENSIONAL ANALYSIS 


OF COUPLED-CAVITY TRAVELING WAVE TUBES 

Thomas A. O'Malley 

Analex Corporation 
Cleveland, Ohio 44135 


SUMMARY 

The use of the coupled-cavity traveling wave tube (TWT) for space com- 
munications has led to an Increased Interest In the efficiency of these de- 
vices. Efficiency Improvements have generally been along two lines: (1) the 

development of efficient multistage depressed collectors to recover a high 
percentage of the kinetic power In the spent beam; and (2) Improvements In the 
efficiency of the basic Interaction process through velocity resynchronization 
and other methods. To pursue the second line, we have recently developed a 
flexible, three-dimensional, axlally-symmetrlc, large-signal computer program 
written In FORTRAN IV for use on the IBM 370 time-sharing system. The present 
report Is a users' manual for this program. 

The basic programming approach consists of dividing the beam Into rings 
and calculating the trajectories of these rings as they pass through a sequence 
of cavities each of which may have different geometrical and electrical proper- 
ties. Extreme flexibility Is provided In the variety of tube features that can 
be modeled since each cavity has Individually entered Input parameters. The 
program can handle lumped or distributed severs, Input and output couplers, 
cavity match details provided for cavities near the end of a stack, voltage- 
jump velocity resynchronization with an arbitrary number of discrete steps (up 
to one per cavity), and velocity taper designs of almost arbitrary complexity. 
Backward waves can be handled by an Interatlve procedure. 


INTRODUCTION 

The use of the coupled-cavity traveling wave tube (TWT) for space com- 
munications has led to an Increased Interest In the efficiency of these de- 
vices. Efficiency Improvements have generally been along two lines: (1) the 

development of efficient multistage depressed collectors to recover a high 
percentage of the kinetic power In the spent beam; and (2) Improvements In the 
efficiency of the basic Interaction process through velocity resynchronization 
and other methods. The NASA Lewis Research Center has been Interested In the 
first approach to efficiency enhancement for sometime (refs. 1 to 3). In the 
second approach, we have developed a flexible, one-dimensional, large-signal 
computer program (refs. 4 to 6). This computer program, because of the limita- 
tions of one-dimensional motion. Is not highly accurate In Its calculation of 
electron trajectories or beam-wave Interaction. Also, one-dimensional motion 
does not allow calculation of beam Interception or vector velocities of elec- 
trons. Beam Interception data Is required In order to determine the properties 
required of an effective magnetic focusing system. Vector velocity data Is 
particularly Important In the design of efficient collectors. The vector 


velocities at the output of the tube are required as Input data to a computer 
program which analyzes collectors. 

Because of the need to overcome the aforementioned limitations Inherent 
In a one-dimensional program, we have recently developed a three-dimensional 
program written In FORTRAN IV. This program retains all of the features In 
the one-dimensional program and Is essentially a generalization of that program 
to a three-dimensional, axial ly-symmetr 1c simulation. 

The basic programming approach Is to follow a ring-model beam through a 
sequence of cavities each of which may have different geometrical and elec- 
trical properties. Extreme flexibility Is provided In the variety of tube 
features that can be modeled since each cavity has Individually entered Input 
parameters. The program can handle lumped or distributed severs. Input and 
output couplers, cavity match details provided for cavities near the end of a 
stack, voltage-jump velocity resynchronization with an arbitrary number of 
discrete steps (up to one per cavity) and velocity taper designs of almost 
arbitrary complexity. Backward waves can be handled by an Iterative procedure 
similar to that described by Kino, et al. (ref. 7). 

In this report, we first discuss the computer model of the coupled-cavity 
TWT. A detailed discussion of the program structure follows. The Input data, 
output data, program variables, and program subroutines are then described. 


MODEL OF COUPLED-CAVITY TWT 

The model of a coupled-cavity TWT Is Illustrated In figure 1. The beam 
Is divided Into a series of disks, and each disk Is subdivided Into an arbi- 
trary number of rings. A complete description of the beam trajectory Is ob- 
tained by following the rings contained In a single beam wavelength (ref. 8). 
The rings may penetrate each other In both the radial and axial directions. 

The rings expand or contract radially according to the radial forces acting on 
them. The rings do not contract or expand In the axial direction; the axial 
thickness of a ring Is assumed to be constant. Also the rings do not become 
warped; that Is, the trace of a ring in the r-z plane Is always a rectangle. 

The axial and radial thickness of a ring comes Into play only In the 
modeling of beam Interception and In the calculation of space charge forces. 
For all other forces, It Is assumed that the charge and mass of a ring are 
concentrated at the ring's centroid radius and that the ring has zero axial 
thickness. 

Each ring, as It enters the tube, has a charge equal to 1/R of the 
charge of the disk, where R Is the number of rings per disk. Also the Ini- 
tial cross-sectional area Is the same for all rings. For the case of R = 3, 
the ring dimensions at the tube entrance are given In figure 2. The axial 
thickness a d of all rings Is X e /N d , where x e Is the beam wavelength and 
N d Is the number of disks In a beam wavelength. 

For an arbitrary value of R, we have 


j = V ^2R b ’ ^ = 1 
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In modeling beam Interception, we must make assumptions on the values of 
the Inner and outer radii, given the value of the centroid radius. The assump- 
tion Is that the ratio of the Inner radius to centroid radius and the ratio of 
the outer radius to centroid radius remain constant. These ratios are deter- 
mined from equations (2) and (3). At the end of an Integration step. If the 
outer radius of a ring exceeds the tunnel radius a, then the outer radius Is 
redefined to be a, and the portion of the ring located beyond r = a Is 
assumed to be Intercepted. The remaining partial ring, since Its outer radius 
has been redefined, has a redefined value of centroid radius given by 



j = 1, ..., R 


(4) 


The program has an option allowing a one-dimensional beam to be used In 
the first N cavities, where N Is an Input parameter, and a switch Is then 
made to a three-dimensional simulation. The results for the first N cavities 
would be nearly Identical to those of the one-dimensional program described In 
reference 4 because the same model Is used In both programs. This report will 
discuss only the three-dimensional model since the one-dimensional model Is 
thoroughly documented In references 4 and 5. A basic understanding of refer- 
ences 4 and 5 Is presupposed In this report. The tube body Is treated as a 
conducting tunnel of radius a divided axially Into a series of discrete 
cells, where the length of the k** 1 cell Is denoted by L^. In the center of 
each cell Is a gap of length 28.^. Impressed across the k*h gap Is a com- 
plex voltage V^e' 10 *. 

In the absence of a beam, there will be a simple known relation among the 
various Vu. For a forward-propagating wave In the pass band of a uniform 
structure, 


-(a+1B. ) L. 

V. v k < 5 > 

where a and Bi are known from cold-test measurements. The power flowing 
along the structure will be given by 


P 


k " 



where also Is known from cold-test measurements. 


( 6 ) 
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The voltage wave In the presence of the beam is calculated by adding a 
complex Induced voltage to the propagating wave, that Is, 



(c+IB, )L. 
e “v**. 


where AVj< + -| Is an Induced voltage to be defined later. 


(7) 


The body of a coupled-cavity TWT Is supposedly a periodic structure. If 
this were the case, all cavities would have the same L, a, oL, B-|L, and Z. 
However, a real tube has Input and output couplers, severs, and perhaps veloc- 
ity tapers and cavity match details, all of which require Individual non- 
periodic treatment In the model. Thus, the properties of Individual cavities 
are separately specified Input parameters In the program, and a great variety 
of tube design variations such as complicated velocity tapers or voltage-jump 
configurations can be easily modeled. Equation (7) Is easily modified to 
handle such a quasl-perlodlc structure. 
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The factor "^Tyy\ /Z^ ensures that power flow Is conserved. 


PROGRAM STRUCTURE 

The program consists of two major elements: the dynamics of the beam 

rings and the beam-wave coupling analysis. These two elements are described 
In detail In the following sections. 


Dynamics of Beam Rings 

Equations of motion . - As discussed In reference 4, page 4, It Is more 
convenient to let the Independent variable be the axial position z rather 
than the time t, and we calculate the functions r(z) and t(z). Here r(z) 

Is the centroid radius of the ring In question when It reaches the axial posi- 
tion z, and t(z) Is the time of arrival of the ring at the axial position z. 
The transformation from equations of motion with t as Independent variable 
to equations of motion with z as Independent variable Is given by 



(8) 
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(9) 
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d 2 z dr 
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Instead of using t, z, and r, we use the normalized variables e, e, and 
p defined by 


9 = wt 

Z 

a 
r 


5 = 


( 10 ) 

( 11 ) 

( 12 ) 


In terms of the normalized variables, the transformation equations are 




( 13 ) 


dfjo _ d£i dp 
d 2 P _ de 2 de 2 
dE 2 = /dl \ 2 


( 14 ) 


Although the normalized variables are used In the program, the unnormallzed 
variables will be used In any explanation for which greater clarity can be 
achieved. At each Integration step, the usual equations of motion with t as 
Independent variable are used to calculate d 2 z/dt 2 and d 2 r/dt 2 . Knowing 
these quantities, we then calculate d 2 t/dz 2 and d 2 r/dz 2 from equations (8) 

-k m A / O \ ' ' ' 


The equations of motion with t as Independent variable may be In either 
relativistic or nonrelatl vlstlc form. The equations for d 2 z/dt 2 and d 2 r/dt 2 
for the ring In question are 


f = B r (E , _ + E - rtJB ) + U i 
dt 2 (m rf,z sc , z r' zj 


( 15 ) 
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rf.r + E sc,r + r * B z> + B r0 ♦ U r 


(16) 


where to Is the angular velocity of the ring and 8 


Is defined as 
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13 = 1 (nonrelatlvlstlc case) 

8 = \j 1 - ^ (relativistic case) 


( 17 ) 


In the above equations, B z and B r are the axial and radial components of the 
magnetic Induction at the location of the ring, E r f z and E r f r are the com- 
ponents of the rf electric field, and E sc>z and £ sc r are tfie components of 
the space-charge field. For the nonrelatlvlstlc case,*U z and U r are zero. 
For the relativistic case, U z and U r are complex expressions Involving a 
factor W defined as 


U = 


2„2 3 jdt (E rf,r + E sc,r } + dt (E rf,z + E sc,z ) 
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where Is the azimuthal magnetic Induction given by 

-vj. .... r 
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max 


In the above expression, p 0 Is the permeability of free space, r max Is the 
outer radius of the outermost ring, and I 0 ,ave I s average beam current. 
The angular velocity to of the ring Is given by 


to = to - * 
o m 


[ 'J'(r.z) - *( r 0 ,z 0 )-| 

^ J 


(22) 


where ^»( r , z) Is the magnetic flux at the ring location (r,z), and the sub- 
script o denotes conditions at the tube entrance. 

The coordinate system used In the program has Its origin at the beginning 
of the first cavity. The k** 1 cavity Is divided Into N z equal parts of 
length az^, given by 


Az k = N z 


The n*h node In the first cavity, denozed by z n -| , Is defined by 
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(23) 


. ( n - l) 4Z, 


n = 1 N. 


The nth node 1n the k th cavity, denoted by z nk , Is given by 


= m E l Lfn + 


("-*) 


T AZ, 


n = 1, 


2 N, 


where N c Is the number of cavities. A point that lies midway between two 

ah' an f?’ k H 3 ? d z n+l,k Is denoted by z n+1/2 k . Such a point Is called the 

k l f1g * 3) * The numerical Integration of equations 
(8) and (9) for the 1 ™ r1ng 1s such that t (z) and (z) and th J f1 ' 

V6S ° f th6S ? f . Unct1ons are evaluated at antinodes! and the second * 

same at f nodes * J he numer1cal Integration scheme Is the 

tW 7 l Thi t iH!n+f«? d ^ refarence equations (8) to (11), for calculating 
tlCZ) • The Identical scheme Is used for calculating r^(z). 

f , numerical Integration method used here Is not as accurate, for a 

ll+hrtH SteP TK 1Ze ^ as . the w1del V used Runge-Kutta methods or predictor-corrector 
!r!!i dS, H The advantage of the method used here Is that It requires evaluating 

aul re^puainaf ? nly 0 " C ? pe [, 1ntegrat1on step, whereas the other methods re- 
qulres evaluating accelerations at least twice per Integration step. To be 

competitive with the method used here, another method requiring evaluation of 

HnMhi er Jh l0n ! tWl ^ 6 Per Step WOuld have to ach1eve comparable accuracy with 
Rnnal^Mtf Ste Jh S H Ze ' ? U ? h 3 com P ar1son between the method used here and a 
f p JJ e ' K “ tta method requiring evaluation of accelerations twice per step has 
been made. The method used here achieved the greater accuracy. 

. ~ At times It Is convenient to use a single subscript 

nf *?th b ? U ?J er; examp1e ’ would denote the time of arrival 

At JihprVim Pl ?+ ? t the aX al pos1t1on z * where 1 runs from 1 to RN d . 

At other times It Is convenient to use a double subscript to denote therlna 

uhi^n e ? ample ’ ^^ z ^th 0u1d denote the time of arrival of the ring 

r^s h f^i 9 i n tn 1y p WaS A Ue < Jth r1n9 mak1ng up the lth d1sk » where J 9 
runs from 1 to R and 1 runs from 1 to N d . J 

The Initial total velocity of all rings Is 
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(26) 
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where h» c Is the cathode flux at an outermost ring and ri>« 0 Is the Initial 
radius of an outermost ring. The Initial axial velocity Is ' 


' J 1.0 = ^ 1.0 


}J = 1 R 

1 = 1 N, 


(27) 


where rjj f0 Is the Initial radius of the (j,1) th ring. From equation (27), 


Zjj Is not the same for all rings. It Is the largest for the Innermost rings 
and the smallest for the outermost rings. It follows that the beam wavelength. 


defined as z.^ divided by frequency. Is not the same for all rings. This 

Is an Inconvenience we wish to avoid, so we assign to z„ q the value It 
would have If R were equal to 1. 


z 


Jl.o 




where b Is the Initial beam radius. 


(28) 


Since tj^(z) and rj^(z) are the trajectory functions to be calculated, 
we must specify the Initial values of these functions and their first deriva- 
tives. These are given by 


O - 1)1, 




U J1 
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dz 
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j = 1. 
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(29) # 

(30) 

(31) 

(32) 


For some applications dm/dz Is nonzero. For such cases, dr<«/dz Is speci- 
fied In the Input data. J1 


Space charge forces . - The starting point In developing the space-charge 
model Is an expression for the potential field due to a point charge q In 
a cylindrical tunnel. From reference 5, equation (10), the potential at 
( r 2»**2* z ) 1s 
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, M >l ms r V , M tl ms r l 1 
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(33) 


where a Is the tunnel radius, e 0 Is the permittivity of free space, and 


J s (p ms a) = °* 


so 


11 for s = 0 
)o for s * 0 


(34) 


and the charge Is located at ( r-j ,0-j , 0) . We can obtain the electric field set 
up by the point charge from the equation E = -W. By a simple Integration, 
we can then find the electric field due to a charged ring of charge q^ and 
radius r^ . The field at (r 2 ,z) Is 


E z (r ? ,z) = 


sgn(z)q ] 

2irc a 2 
o 



J o (tJ m r l )J o (ti m r 2 ) ^m' 21 




(35) 
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r 



.z) = 



J (p r.)J.(p r_) ~* J m ,z * 

o' H m 1 ' l' H m 2' e 




(36) 


where p mo has been replaced by p m . 

The space charge force on a reference ring of charge q 2 and radius r 2 
and located at z, due to a source ring of charge q-| and radius r-| and 
located at z = 0, Is then given by 

f 21,z (r l -r 2’ Z) = q 2 E z (r 2 ,Z> (37> 

r 21.r <r r r 2, z) ' q 2 E r (r 2’ Z) (38 > 

Equations (35) to (38) cannot be used to calculate space-charge forces 
because the summations will approach Infinity If r^ approaches r 2 and z 
approaches zero. We must take Into account the nonzero axial and radial thick- 
ness of the rings. A fully general calculation can be made by Integrating over 
both the source ring, located at (r^.O), and the reference ring, located at 
(r 2 ,z). Because of the complexity Involved, It Is not feasible to calculate 
these expressions each time the space charge force on one ring due to another 
Is desired. Instead the radial and axial space-charge forces are calculated 
for an array of values of r^ , r 2 , and z and stored In a table. This cal- 
culation Is done once, just prior to the beginning of the simulation. In the 
simulation, the space-charge forces are evaluated by a three-dimensional linear 
Interpolation on the table. 
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In carrying out the Integrations over the axial and radial thicknesses of 
the rings, assumptions must be made on what these thicknesses are. It Is 
assumed that the axial thickness, st d = x e /N d (see fig. 2), remains constant. 

The radial thickness Is determined from the assumption that the ratio of the 
Inner radius to centroid radius and the ratio of the outer radius to centroid 
radius remain constant. The Inner and outer radii are obtained from a given 
centroid radius by using equations (2) and (3). Since these radii depend on 
the ring type j, the space-charge force Is really a function of five vari- 
ables: r-| , r 2 » z, J] , J2» Since j] and j2 can take on all values from 1 to 

R, R^ tables must be calculated. The storage and time requirements for cal- 

culating this many tables may be excessive for some applications. These re- 
quirements are reduced If a less general formulation Is adopted. We note that 
when calculating the rf electric-field forces on a reference ring, the approach 
we take Is to calculate this field and evaluate It at the centroid values of r 

and z of the reference ring. We can adopt the same approach with the space- 

charge field forces; that Is, we calculate the space charge field by Inte- 
grating over the axial and radial thicknesses of all source rings and then 
evaluate the field at the centroid values of r and z of the reference ring. 
Since we Integrate over source rings but not reference rings, the space-charge 
force Is now a function of four variables: r-i , r 2 , z, ji. In this approach 

only R tables must be calculated. 


The two models for calculating space-charge forces described In the previ- 
ous paragraph will be referred to as the R 2 -tables model and the R-tables 
model. We will discuss the relative merits of the two models. But first we 
mention a modification that can be made to the R-tables model. The modifica- 
tion Is to Integrate over the axial thickness of the reference ring. This 
modification only slightly Increases the computational complexity and does not 
Increase the number of required tables. Also space-charge forces become 
smoother functions of z. To Indicate the Increased smoothness as a function 
of z, figure 4 shows the axial space-charge force on an Innermost ring due to 
another Innermost ring, as a function of z, with and without the modification 
just described. Another advantage of this modification Is that, for cases 
where n = 0, r 2 > 0 and z < a d , the summation for the radial space-charge 
force Is a convergent series, whereas this series Is divergent without the 
modification. Also, the summations for both axial and radial forces converge 
more rapidly with the modification. In all further discussions of the R-tables 
model, It Is assumed that the modification Is made. 

Comparison of the R 2 -tables model to the R-tables model for a variety 
of cases shows that disagreement In the space-charge force on a reference ring 
due to any one source ring Is less than 7 percent In practical situations. 
Typically there Is considerably less disagreement In the space-charge force on 
a reference ring due to al_l_ source rings. Figure 5 shows a typical comparison 
of the R -tables model to the R-tables model. We consider these disagreements 
to be reasonably small, especially when compared to the error already present 
In the model from the assumption that there are no gaps In the tunnel. We 
therefore adopt the R-tables model because of Its reduced computational and 
storage requirements. 


The expressions for the axial and radial space-charge forces for the 
R-tables model are obtained by carrying out the appropriate Integrations In- 
volved In equations (35) and (36). Letting r,., n and r, out be the Inner 
and outer radii of the source ring, we obtain 
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where the functions G a , Gj,, G c , G<j are defined as 
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For the region r-| = 0, r 2 > 0, z < 8,^, the summation for the radial 
space-charge force Is a very slowly converging series. For this region, we 
use an alternate computation. We subdivide the region Into two regions: 

(1) r^ = 0, 0 < rg < 0.125a, z < 1^* and (2) r l = r 2 > 0.125a, z < 
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For the first region the tunnel wall Is far removed from the source ring and 
reference ring, and we consider the rings to be In free space. For this case 
the radial force on a reference ring of zero axial thickness Is given by 




The radial force on a reference ring of axial thickness 8. d Is obtained by 
numerically Integrating the above equation. For the second region, the pres- 
ence of the tunnel wall cannot be neglected. It can be shown that 



Thus F21,r( r l = O.^.z) should differ only slightly from F 21 r ( r l = Ar sc» r 2* z )* 
where r-| = Ar sc represents the first nonzero r-| grid point! We use this 
approximation In the second region for calculating F 2 i tr (rj = 0 ,r 2 ,z). 

F 21,r (r l = °’ r 2’ z) = F 21,r (r l = Ar sc ,r 2 ,z) ’ r 2 > °- 125a (43) 


For the region r-j = 0, r 2 = 0, z < a, d , the summation for the axial space- 
charge force Is a divergent series. For this case, the source and reference 
"rings" are actually line charges of length 8 d which partially overlap one 
another. We choose to truncate this summation after 100 terms, a somewhat 
arbitrary choice, and use the result for the axial force. It should be noted 
that the grid points In this region will be used rarely. If at all. They will 
be used only when a reference ring and a source ring, at the same time, have 

centroid radii less than Ar sc and are separated axially by less than a, d . 

From figure 4, we have an Idea of the smoothness of the axial space-charge 

force as a function of z. The radial force as a function of z has a similar 

degree of smoothness. From examining these functions we can determine the 
spacing between grid points In the z-dlrectlon that Is required In the space- 
charge tables for sufficiently accurate Interpolation. In figure 4, a d /a = 1/3. 
The shape of these functions varies somewhat with varying values of a, d /a. 

After examining these functions for a range of values of i d /a, we decided to 
use two different grid spaclngs. Letting az sc be the spacing between z grid 
points, we use 
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(44) 


A Z 



d’ 


0 < z < 38. . 
- - d 


A Z 



z > 38. 


(45) 


Let z max he the last z grid point In the space-charge tables. In the 
simulation, If the axial separation distance between two rings exceeds z max , 
the space-charge forces are set to zero. We calculate z max such that the 
error In setting these forces to zero Is less than 1 percent. An empirical 
formula Is used to obtain z max . 


In 


max 


= 0.09 


0.345 1 

,V a ) 


In 


0.693 


- 0.26 


' 0.345 1 


0.693 


+ 2.41, 


-f- > 0.127 (46) 

d 


Z 

max 

a 


2 . 22 , 


— < 0.127 
a 


(47) 


It Is undesirable for z max to be nearly equal to an Integer multiple of a.^. 

To avoid such an occurrence, z max Is Increased slightly so that z/a.^ Is equal 
to an Integer plus 1/2. To cut down on storage requirements, we use no more 
than 26 z grid points. If Zmax^d greater than 9 1/2, we must Increase 
the grid spacing Az sc to satisfy the constraint of a maximum of 26 z grid 
points. When z ma x /a d ls greater than 9 1/2, Az sc Is redefined to be 


A z 



0 < z < 31 . 
— — d 


(48) 


AZ 


SC 


z 

max 
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Z > 3 ft. . 
d 


(49) 


We should also examine the space-charge forces as functions of either of 
the two r variables. We concentrate our attention on the radial force since 
this Is the more rapidly-varying force. We consider an example of the func- 
tional dependence on r In figure 6. For this case, R = 3 and the source 
ring Is a "middle" ring of Inner radius 0.326a, centroid radius 0.4a, and outer 
radius 0.462a. Figure 6 Is the radial space-charge force on a reference ring, 
due to the source ring, as a function of the reference ring's centroid radius. 
The two rings are at the same axial position. The vertical dotted lines re- 
present the Inner and outer radius of the source ring. The other dotted line 
represents the linear Interpolation of this function for a radial grid spacing 
of Ar sc = a/N sc>r with N sc r = 10. The accuracy of the Interpolation Is 
poor In the regions near r =’0.326a and r = 0.462a. At r = 0.462a the 
Interpolation error Is 27 percent. When N sc r Is Increased to 20, this error 
Is reduced to 12 percent. For rough-cut runs! N sc r = 10 may provide suffi- 
cient accuracy, but N SCjr = 20 should be used when greater accuracy Is de- 
sired. The larger value’of N sc>r should certainly be used In cases where the 
average beam diameter Is small compared to the tunnel diameter. 
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The radial space charge force becomes much more rapidly varying as either 
of the r variables becomes smaller. For this reason, we use a smaller grid 
spacing for small values of r In such a way that the total number of grid 
points Is still N sc r + 1 . For larger values of r, the spacing between grid 
points Is only slightly greater than It would be for completely uniform spacing. 
The radial grid spacing Is given by (we assume N sc>r Is a multiple of 5) 


Ar 


sc 


5a 


sc,r 


0 - r - iJ 


(50) 


sc,r 


Ar 


5a <"sc.r -’> 


SC 


r > 


4N 


(51) 


sc,r 


sc,r 


The space charge force on a given ring may be a slowly varying function 
of the Independent variable z. This Is likely to be true In the beginning 
cavities of the tube. To take advantage of the slowly varying nature of the 
space charge force, the calculation Is not done on every Integration step. 

The calculation of space charge forces Is Initially done every step, 

where N Is an Input. The value of N may change during the simulation In a 
manner to be described later. On those Integration steps where space charge 
forces are not calculated, they are approximated by quadratic curve fits. 

Each ring has Its own curve fit. These curve fits are obtained by fitting a 
quadratic polynomial through the latest three calculations of space charge 
forces. Let us say that these three calculations were done when z was equal 
to z a , Z|j, and z c . At an Integration step where a new space charge force 
calculation Is due, at z = z^, the values calculated are compared with those 
obtained when the curve fits based on the points z a , z and z c are eval- 
uated at Z(j. If the difference In the values for any ring Is exceeded by a 
specified tolerance, the next calculation of space charge forces occurs 1/2 N 
steps later. In general, N Is halved each time the specified tolerance Is 
exceeded. 


For calculating space charge forces, the positions of all rings must be 
known at a time t. This Information Is not readily available since z rather 
than t Is the Independent variable. To calculate the space charge force on 
the 1 th ring located at some z, we must determine the positions of all the 
other rings at the time t^(z). Consider the j th ring. The time of arrival 
tj(z) at z Is known. By the periodicity of the motion, the ring that origi- 
nally was m wavelengths behind the j th ring arrives at z at the time 
tj(z) + mT, where T Is the reciprocal of frequency. Similarly, the ring 
that originally was m wavelengths ahead of the ring arrives at z at 

the time tj(z) + (-m)T. We consider all such rings and find m such that 

|t^(z) - (tj(z) + mT) | = Minimum (52) 

We may call the ring In question the ring In the m th cycle. We assume 

that the j*" ring In any other cycle Is too far away from the 1*h ring, at 
the time t^(z), to contribute significantly to the space charge force on the 

1 th ring. Let (r^z™) be the position of the ring In the m th cycle 

at the time t^(z). We obtain ( , zT) by assuming that this ring has con- 
stant axial velocity vj z (z) and constant radial velocity vj r (z) In the time 
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Interval spanned by t^(z) and tj(z) + mT. Then (rJ,Zj) Is given by 

r J = + v jr^ z ^ t 1^ z ^ " + mT ^ (53) 

Z j = Z + v jz (z)[t 1 (z > ~ *V Z * + mT ^ (54) 

The space charge force Is then calculated based on the reference ring radius 
ri(z), the source ring radius rj, and the axial separation distance z - Zj. 
The calculation Is done for all rings, excluding the 1*h ring, and the re- 
sults are summed to obtain the total space charge force on the 1*h ring. 

When the axial separation distance z - z? Is greater than z^x (see eqs. 

(46) and (47)), the space charge calculation Is not done for this term In the 
summation. 


As discussed In reference 2, relativistic corrections are applied to 
axial lengths. In the space-charge equations, z and 8.^ are replaced by 
z* and 8^ given by 


8 


d 


z' 



(55) 


(56) 


Radiofrequency electric field forces . - If the complex voltage V on a 
given gap Is known, the electric field In the neighboring region can be written 
In an expansion of the form 


where 


F (r 7 t) - -- ^ 

z' ’ ’ > 28 slnh(p) LI ( Ym a) 
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(60) 


a 


ID 



+ 
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(61) 


and 0 O L Is the lowest order phase shift per cavity. The shaping factor y 
Is an Input parameter. 

We define Q z (r,z) and Q r (r,z) as 


Q z (r,z) 


slnh(y) 



LI o<V> 



(62) 


CO 

q r (r,z) ’ slnh(v) 


IVmVV) -'V 

’ r m U o ( V> 


m = -oo 

Equations (57) and (58) can now be written as 

E z (r.z,t) = Q z (r,z) e 1wt 


E r (r,z,t) = Q r (r,z) e 1wt 


(63) 


(64) 

(65) 


Equations (64) and (65) are appropriate to a forward-traveling circuit wave. 
For a combination of a forward- and backward-traveling wave, the appropriate 
expressions are 


E z (r,z,t) 


= [Q z (r,z)V f + Q*(r ,z)V b ]e^ 0,1 


( 66 ) 


E r (r,z,t) = ^ [Q r (r,z)V f + Q*(r,z)V b ]e 1wt (67) 

The electric-field force on a ring of charge q and located at (r,z) at time 
t Is then given by 

F Z = R e [fl «M r ’ z)V f + Q*(r,z))e Wt ] (68) 

F r = R e[f? <V r ’ z)V f + Qf(r.z))e Wt ] (69) 

Many of the quantities In equations (57) to (69) change from cavity to cavity. 
When a particular cavity, let us say the k th cavity, Is being discussed, the 
subscript k will be appended to these quantities. We now consider the I** 1 
ring In the k th cavity located at the n th node z^. The centroid radius 
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of the ring Is r^(z n |<). The electric-field force Is evaluated at r = n(z n |<) 
and z = z^. From equations (68) and (69), the electric-field force on the 
ring Is 


F z,1nk “ R e 


F r,1nk R e 


\_2l k (Q zk (r 1 ( 
[ 21 ^ (Q rk (r 1 


< z nk>’ z nk> V fk + ^zk^ r 1^ z nk^ ,z nk^ V bk^ e 




(z nk ),z nk )V fk + Q rk (r 1 (z nk ) ,z nk )V bk )e 


'] 

1< z nk>J 


(70) 


(71) 


The Q Z | c (r,z) and Q r | C (r,z) functions are evaluated for an array of r and 
z values and stored In tables when cavity k Is entered. The z values are 
the nodes z n |<, n = 1, N z . The r values are equally spaced values 
from r = 0 to r = a. The number of r values Is an Input parameter. The 
Qzk( r l( z nk) * z nk) and Q r k( r l( z nk) » z nk) equations (70) and (71) are evaluated 
by doing a linear Interpolation In r on these stored tables. The number of 
terms Included In the summations for Q z |< and Q,^ Is determined by an Input 
parameter. If cavity k has the same L^, 8,^, 8^, p^ as cavity k-1 , then 
Q z , k-1 = Qzk anc * Qr,k-1 = Qrk anc * new tables do not have to be calculated. 
When new tables do have to be calculated, the old tables are no longer needed, 
and the new tables may occupy the same storage space as the old tables. The 
electric-field forces are evaluated by equations (70) and (71). 

Voltage-lump electric field forces . - If a dc voltage Vj Is Imposed 
across the gap of a cavity, the resultant electric field In the cavity Is given 
by (ref. 5) 


where 


m 


,(r,z) = -f 1 * JC, Vo< k m r)cos( V>] 

- xfjf, Vi<v> s1n <y)] 


2p[p slnh p cos k 8. + k 8. cosh p sin k a] 
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We define R z (r,z) and R r (r,z) as 


(72) 


(73) 


(74) 


(75) 


R z < r .*> - 1 * e wy cos <v> 

m = 1 


(76) 
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( 77 ) 


R r (r - Z) ’ £, WV> s1n < k m z > 

ID = I 

Equations (72) and (73) can now be written as 


* ' f 

E z (r.z)= 

(78) 

-V n R,(r,z) 

E r <r.z) . ' 

(79) 


As before, when referring to the cavity, the subscript k will be ap- 

pended to those quantities which may change from cavity to cavity. We now con- 
sider the 1™ ring In the k^ cavity located at the n*h node z n |<. 

From equations (78) and (79), the electric-field force on the ring Is 


~ q 1 V Jk R zk (r 1 (z nk ),z nk ) 


z,1nk 


(80) 


r,1nk 




(81) 


The R Z | c (r,z) and R r k(r,z) functions are calculated and stored In a table and 
have the same features as the Q Z |<(r,z) and Q r j<(r,z) for rf electric-field 
forces 

Magnetic field forces . - The program provides options for simulating two 
types of magnetic focusing fields: (1) uniform solenoldal focusing and (2) 

periodic permanent magnetic (PPM) focusing. 

The magnetic field for uniform solenoldal focusing Is given simply by 


B z ( r , z) = B 0 

(82) 

B r ( r , z) = 0 

(83) 


where the constant B 0 Is a program Input. 

The second type of magnetic focusing Is PPM focusing. We consider first 
"single period" PPM focusing. Figure 7 shows the geometry of the ring magnets 
which are centered along the z-axls. The gap length between magnets Is 2g, the 
Inner diameter Is 2a, and the magnetic period is P. Also shown In figure 7 Is 
the shape of the axial magnetic field at r = a. The origin Is at the center 
of a gap so that B z Is an even function of z. The assumed shape of B z (a,z) 
In the region |z| < P/2 Is 


B z (a,z) = B q cosh yz, 0 < |z| < | 


(84) 
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(85) 


B z (a,z) = 0, 


f < | 2 | < j - g 


B z (a,z) = -B Q cosh 


[=H)]- 


I- s < III <f 


( 86 ) 


where the shaping factor p Is an Input parameter. 


The magnetic field can be found using a Fourier series expansion. We note 
that B z Is axially symmetric and satisfies LaPlace's equation, and B r Is 
zero at z = ± P/2. The Fourier series expansion which satisfies these condi- 
tions and has the required behavior at r = 0 Is 


B z (r - Z) ■ n £, a n I o (k n r) tos(k r, z> 


(87) 


where k n Is given by 


2nw 


n “ P 

Using the equation dlv § = 0, we can show that B_ Is given by 


( 88 ) 


V.D - ^ a „W> cos(k n z) 


(89) 


n = 1 


The coefficients a n can be evaluated If B z (r,z) Is a known function of z 
at some value of r, let us say r = r 0 . We obtain 


a_ = 


Pljk T ) 
o n o 


/ 

VCi 


B z<V z > 


cos 


k n z 


dz 


(90) 


If r 0 = a, we use the assumed field shapes given by equations (84) to (86). 
Substituting Into equation (90), we obtain 

8B Q [p slnh pg cos k n g + k n cosh pg sin k R g] 

3,1 p (“ 2 * 3I)W> 

a_ = 0, n even 
n 


In practice, the field shapes at r = a are not precisely known and cannot be 
conveniently measured. Alternatively, the on-axis values B z (0,z) are con- 
veniently measurable and can be used In equation (90), letting r 0 = 0, to 
evaluate a n . For the case of measured data, the Integral In equation (90) 

Is evaluated numerically. 


•\ 


n odd 


(91) 
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"Double period" PPM focusing Is treated In much the same way as "single 
period" PPM focusing. Figure 8 shows the geometry of the ring magnets which 
are centered along the z-axls. As before, the gap length between magnets Is 
2g, the Inner diameter Is 2a, and the magnetic period Is P. Also shown In 
figure 8 Is the shape of the axial magnetic field at r = a. Whereas the 
single periodic case has Its origin at the center of a gap, the double periodic 
case has Its origin midway between two adjacent gaps. The origin Is placed 
here In order to make B z an even function of z. The assumed shape of 
B z (a,z) In the region |z| < P/2 Is 


B z (a,z) =0, 0 < |z| < g - g 


B,(a,z) = B^ cosh 

Z 0 


[4 ■ !)]• 


o-g<UI<o + g 


B z (a,z) =0, I + g < |Z| < - g 


B z (a,z) = -B Q cosh 


[4-r)]- 


B z (a,z) =0, jp + g < |z| < | 


(92) 

(93) 

(94) 

(95) 
(94) 


In the same way as before, the magnetic field can be found using a 
Fourier series expansion. The results for B z (r,z) and_ B r (r,z) are 
Identical to those In equations (87) to (90). If r 0 = a, we use the assumed 
field shapes given by equations (92) to (96). Substituting Into equation 
(90) , we obtain 


a_ = 


(2 cos -p) (8B q ) [p slnh vg cos k p g + k n cosh pg sin k n g] 


* k n) 


n odd 


a n ■ “• 
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(97) 




As before, the on-axis values B z (0,z), obtained from measured data, may be 
used In equation (90), letting r 0 = 0, to evaluate a n . 

The magnetic flux t|»(r,z) for both "single period" and "double period" PPM 
focusing Is given by 


4 »(r,z) 
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The magnetic flux Is needed for calculating the angular velocity 0 of a ring 

• • 

(see eq. (22)). Knowing 0 , the axial magnetic force r0B and the radial 
• 8 

magnetic force r0B z on a ring can be calculated. 

As before, when referring to the k* h cavity, the subscript k will be 
appended to those quantities which may change from cavity to cavity. We now 

consider the 1^ ring In the k*^ 1 cavity located at the n** 1 node z^. 

The angular velocity 0 ^ of the ring Is 


“ink " 


q t ¥ r 1 (z nk>-W - *o 
J ‘ r 1< z nk> 2 


The axial and radial magnetic forces on the ring are 


(99) 


F z,1nk = q 1 r 1 (z nk )0 1nk B rk (r 1 (z nk ),z nk } 


( 100 ) 


r , Ink 


= _q 1 r 1 (z nk )0 1nk B zk (r 1 (z nk ),z nk ) 


( 101 ) 


The B Z | < (r,z), B^r.z), and ^(r.z) functions are calculated and stored In a 
table and have the same features as the Q Z |<(r,z), Q r |<(r,z), R Z |<(r,z), and 
R r k( r . z) tables. 

Rings with small radius . - In this section we discuss the behavior of 
rings which become very small In radius. As we shall see, a modification to 
the simulation will be required. It Is Instructive to analyze the equations 
of motion and obtain an approximation to the radial motion of such rings. 


We look at the case of an Innermost ring, which Is really a disk, since 
this ring type Is the one most likely to become very small. We use non- 
relatlvlstlc equations and assume the cathode flux Is zero. From equation 
(16), the radial equation of motion Is 


0 r 
dt 2 
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c + E 
rf,r sc,r 


I VO J 


( 102 ) 


From equation (22), 0 Is 


0 = za * (r,z) 


m 


2irr 


(103) 


For small r, B z (r,z) can be approximated by B z (0,z). Then «|»(r,z) can be 
approximated by 

*(r,z) = irr2B z (0,z) (104) 
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Substituting equation (104) Into equation (103), we obtain 


-qB (0,z) 

V- 


(105) 


Substitution of equation (105) Into equation (102) yields 


d ? r -Q^<0-*>r 

dt 2 am 2 + m rf.r + E sc,r ) 


(106) 


We transform the above equation to one with z as Independent variable, as- 
suming that the axial velocity v z Is constant over the region of Interest. 


'r -<n^(0,z)r 


4m 2 v 2 


2 (E rf,r + E sc,r ) 


(107) 


For r small enough, the space-charge field Is dominated by the self-force. 
We assume r Is small enough that the space-charge field can be approximated 
by the self-force only. From Gauss' law, the self-force Is closely approxi- 
mated by 


sc,r 2 ire Jt.r 
o d 


(108) 


Equation (107) becomes 


:?!l (0 ’ z)r ^ c , / o 2 

4m V z mv 2 rf ’ r U^mv 2 ^ 


(109) 


The first two terms on the right hand side of equation (109) are negligible com- 
pared to the third term for small enough r. We then obtain the approximation 


dr K 
j_ 2 = r 


(HO) 


K 



2ire mv 8. , 
o z d 


I q/m 1 1 


4ire U V 
O 0 z 


( 111 ) 


Equation (110) cannot be solved analytically; however, we can determine 
the minimum value of r. For shorter notation, we use a "prime" to denote 
differentiation with respect to r. We consider the situation where r(z) and 
r'(z) are known at some point z , . 
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r ( z 1> ’ r l 


r ■ ( z, ) 


r 
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We assume further that r(z) Is decreasing and reaches a minimum value r 2 at 
some point zg. The following derivation for r 2 Is self-explanatory. 


r" 
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2r ' r" 


2Kr 1 
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d /drV 2 K dr 

dz \dz j ~ r dz 



<r') 2 - (rj) 2 


2K ln f7 


Continuing the derivation and noting that r£ = 0, we have 




-(r') 2 /2K r^ 
-(rj) 2 /2K 


We consider a numerical example from a computer run simulating the Com- 
munications Technology Satellite (CTS) 200-watt TWT (ref. 9). The pertinent 
data are: 


r 1 = 0.05a = 0.3175xl0' 4 m 
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(corresponds to 3.5° angle) 


= v = 0.63X10 8 — 
z sec 

I =0.07 amps 
o 


r 2 -3 

— ■ = 0.73x10 J 

d 

The actual trajectory, obtained from numerical Integration of equation (110), 

Is shown In figure 9. In the neighborhood of the minimum point, r" Is changing 
very rapidly and very small Integration steps need to be taken. Since the 
program Is designed to have a fixed number of equally-sized steps per cavity. 

It Is not possible to accurately calculate the type of trajectory shown In 
figure 9, without using a prohibitively large number of steps per cavity. We 
notice from figure 9 that the ring approaches the z-axls so closely that It 
almost appears as though the ring Is reflected from the axis. In the program, 
we simulate such a trajectory as though the ring Is Indeed reflected from the 
z-axls. Assume a ring has radius r and first derivative r' at the end of 
an Integration step and that the step size Is Az. If the ring satisfies the 
following condition. 




With this data, we obtain for 


r + 2Az r* < 0 

then the decision Is made that the normal Integration mode will not accurately 
calculate the ring radius over the next step. Instead the acceleration on the 
ring Is set to zero and the simulation proceeds. When the Integration step Is 
reached at which the ring radius would otherwise become negative, the reflec- 
tion takes place and dr/d z changes sign. After the reflection takes place, 
the actual value of the ring acceleration Is restored, and the normal mode of 
simulation proceeds. In figure 10 the resultant trajectory (dashed line) Is 
compared to the actual trajectory. 

Beam-Wave Coupling Analysis 

The fields Induced by the beam current can be conveniently represented by 
Induced gap voltages. The analysis for obtaining expressions for the Induced 
voltages Is given In reference 5 (pp. 11-15). From equations (66) and (67), 
the forward and backward components of the vector electric field can be 
expressed as 


-V 

£ f (r,z,t) = ^ (r,z) e1wt 

(112) 

-V 

E b (r,z,t) = 2 ^" Q*(r,z) e 1o>t 

(113) 
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where 3(r,z) Is defined as Q z (r,z)a z + Q r (r,z)a r . Equations (46) and (47) of 
reference 5, when generalized to the three-dimensional case, become 

AV f = U J ^*(r,z) • J 1 (r,z)d» (114) 

AV b = fit f ^(r.z) * ^(r.z)dv (115) 

Jv 

where 3i(r,z) Is the fundamental Fourier component of the beam current density 

3(r,z,t), v Is the volume of the cavity In question, and Z Is the cavity's 
Interaction Impedance. 

Since |Q r (r,z)| « |Q z (r,z)| and |J lr (r,z)| « Ui z (r,z)|, we may approximate 
($(r,z) • Ji ( r , z) by Q z (r,z) Ji z (r,z). Appendix B gives a deviation of the 
following expression for J lz (r,z) 



(116) 


In the above expression 
by 


g(r;r. , ,r. .) Is the rectangular function defined 

* 1,1n 1,out 


g( r;r. . ,r, . ) = 1 , 

1 , In 1,out' 
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’ 1,1n 1,out 7 


r, . < r < r. . 

1,1n 1,out 


r < r, , or r > r. , 
1 , 1 n 1,out 


(117) 


where r. , and r, . are the Inner and outer radius of the 1 
1 , In 1 , out 

The backward Induced voltage, from equation (115), now becomes 
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ring. 
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(118) 


Using equation (116), equation (118) becomes 
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(119) 
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The Integration In r can be shown to be 
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We define S(r,z) as follows: 


( 120 ) 
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Equation (120) becomes 

r r \,o^ z) 

* Vln < z > 


O z (r.z) r dr = a r 1f0ut (z) S(r, out (z),z) 


- a 


r 1.1n (z) S < r 1,1„< z >- z > 

With the aid of equation (122), equation (119) becomes 




where Gfc(z) Is given by 


N d R 


: 

1 2 


G b (z)dz 


6h(*) = 77 7 


»E 


G)S,_ 


dd(z) 


sin 


[ r 1.out (z > - r 1,1n< z >] 


2v^(z) -lot^z) 


wit. 


2v.(z) 

1=1 1 
x 2watr 1,out (z) S(r 1,out (z) ’ z) " r l,m (z) S(r 1,1n (z) ’ z) l 


( 121 ) 


( 122 ) 


(123) 


(124) 
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Similarly, AVf Is given by 


AV f - 


where Gf(z) Is given by 


V 


v z > • k T 


E 


1 = 1 



q 1 < z) 


(z) dz 


0)1 „ 


sin 


2v 1 (z) -Iwt^z) 


- r M" U) ] ^ 


2v 1 (z) 


(125) 


x 2 * atr t.out< z > W z,>z) - r 1,ln <2) s *< r 1, t n< 2 >- 2 ” (,26) 

In equations (123) to (126), the Integration variable z Is the axial position 
relative to the center of cavity k. To avoid excessively cumbersome notation, 
the subscript k has been omitted from these equations. The Integrals In 
equations (123) and (125) are evaluated with the use of Simpson's rule with end 
corrections. The end corrections were derived with the assumption that the 
number of nodes per cavity N z Is a multiple of four. Equations (123) and 
(125), with the subscript k restored, now become 


N 


AV 


fk 


* z k Y, wv 


(127) 


n = 1 


AV 


bk 


4z k 4 n G b (z nk> 


(128) 


where, from Simpson’s 

rule with end corrections, 

the 6 n 

are 

. 1 17 

*1 = 3 + 24 

4 n - I- 4 S n i N z 

- 2, If 

n 

Is odd 

x _ 1 2_ 

°2 - 3 24 

4 < n < N, 

- 2, If 

n 

Is even 

i _ 2. + 2— 

*3 - 3 24 

i _ 1 + L. 

d N z -l " 3 24 



1 



(129) 
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From equations (124) and (126), the evaluation of G f (z nk ) and G b (z nk ) re- 
quires the evaluation of S k (r, out (z nk ) z„ k ) and S k (r, )n (z nk ),t nk l for 
1 = 1, ...» NjR. The S k (r,z) function Is calculated and stored In a table In 
the same manner and with the same features as the R zk (r,z), R rk (r,z), Q, k (r,z) 
and Q rk (r, z) tables. The evaluation of S k (r,z) In equations (124) and (126) 
Is done by linearly Interpolating the table In the r variable. 


Calculating the effects of the backward wave requires an Iterative proce- 
dure. In the first pass through the tube, the AV bk are calculated and stored 
for each cavity k. From the AV bk , the backward voltages for a cavity chain 
from cavity kl to cavity k2 are obtained from 


V b,k2 = AV b,k2 


(130) 



AV bk + V 


b,k+l 



-KBi L ) k 


k = k2 - 1, .... kl (131) 


By a cavity chain, we mean a sequence of consecutive cavities whose backward 
voltages are to be calculated. After V bk are known for the desired cavities, 
a second pass Is made through the tube. The second pass yields a new set of 
AV bk that can be used for calculating a set of V bk for a third pass. The 

process continues until convergence Is obtained. In many applications the 

Iterative procedure will converge faster If the Iteration Is done on one cavity 
chain at a time. This would be true, for example, for a tube having severs. 

The program has the capability of performing the Iterative procedure for an 
arbitrary set of cavity chains. 

The computational procedure In the beam-wave Interaction process Is as 
follows: 

(1) When cavity k Is entered, calculate and store the S k (r,z) table 

given by equation (121). If cavity k has the same L k , l k , B ok , p k as 

cavity k— 1 , then S k _-|(r,z) = S k (r,z) and a new table does not have to be 
calculated. 

(2) Obtain a first approximation to AVf k by assuming that the rings 
have constant axial and radial velocities In cavity k. Thus 


n = 1, .... H z (132) 


Knowing r 1 (z nk ), we can readily calculate r iout (z nk ) and r 1 1n (z nk ). 
Knowing t^(z nk ), r^ ou ^(z nk ) and r^ ^ n (z nk ),*we can calculate’the first ap- 
proximation to AV fk from equations (126), (127) and (129). 

(3) Let the forward voltage for cavity k be 


dt(z ) 

V z nk> " V z lk> + -6T~ < z nk - z ,k> 


dz 


dr <z,k> 


r 1 (z nk> ’ r 1< z !k> ♦-dT 1 -< z nk - z lk>. 


dz 


V 


fk = 





+ A V 


fk 
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If there Is a backward voltage, we assume It Is known from equation (131) using 
the set of AVfok obtained from the previous pass through the tube. 

(4) Proceed with numerical Integration of the equations of motion to ob- 
tain trajectory data at each node and antinode In the cavity. 

(5) At the n th node In the cavity, after t^(z n k) and rj(z n k) have 
been obtained, calculate and store 4 n Gf(z n k) and 4 n Gb( z nk)* These quanti- 
ties are the n^ terms In the summations of equations (127) and (128). 

(6) When the last Integration step In cavity k has been done, calculate 
a better approximation to AVfk from equation (127). The AVfk so obtained 
replaces the old AVfk calculated In step 2. Also, calculate AV[>k from 
equation (128). The AV^ are stored for use In the next pass through the 
tube. 

(7) If additional accuracy Is required, make another pass through the 
cavity, repeating steps 3 to 6. 

(8) Repeat steps 1 to 7 for cavity k + 1. 

(9) When the last cavity Is done and the effects of the backward wave 
have to be determined, make a second pass through the tube. The backward volt- 
ages for the second pass are determined from the AV^k of the first pass. 

(10) Hake as many passes through the tube as are required to obtain 
convergence. 


Summary of Program Steps 

The entire computational procedure Is given In flow chart form In figure 
11. The procedure Is summarized by the following steps: 

(1) Read Input data. 

(2) Calculate the tables for space charge forces. 

(3) Begin numerical Integration of equations of motion with z as the 
Independent variable. 

(4) When a new cavity Is entered, let us say the k™ cavity, print 
data for cavity k - 1. If the parameters of cavity k are different from 
those of cavity k - 1, calculate the tables that are required for evaluation 
of rf forces, voltage-jump forces, magnetic forces, and Induced voltages. 

Obtain a first approximation to AVfk by assuming that the rings have con- 
stant velocity throughout cavity k. Attenuate and phase shift Vf k_i and 
vectorial ly add the result to AVfk to obtain the forward voltage for 
cavity k. 

(5) When the last Integration step In cavity k has been done, calculate 
a better approximation to AVfk, an( * use new approximation In place of the 
old. Also, calculate and store AVj,k* 

(6) If additional accuracy Is required, make a second pass through 
cavity k. 

(7) Repeat from step 4 for cavity k + 1. 

(8) If the effects of the backward wave are to be determined, make a 
second pass through the tube. 

(9) Make as many passes through the tube as are required to obtain 
convergence. 


DESCRIPTION OF INPUT DATA 

The following data are required by the program. In the following defini- 
tions, a "magnetic section" Is defined to be one-half of the magnetic period 
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for single-period PPM focusing and Is defined to be one-fourth of the maqnetlc 
period for double-period PPM focusing. 


Name Symbol 


Description 


ACM a Tunnel radius, cm 

ALPHL(K) (oLJfk Voltage attenuation of forward wave In going from cavity K 

to cavity K + 1, dB per cavity; K = 1, LASTCV 

ALPHLR(K) (aL)^ Voltage attenuation of backward wave In going from cavity 

K + 1 to cavity K, dB per cavity; K = 1 , LASTCV. (ALPHLR 
needs to be loaded only If KLOSS = 1.) 

B °( K ) B ok Value of axial magnetic field B z (r,z) at r = a and 

z = midpoint of gap In K^h magnetic section; webers per 

square meter; K = 1, LASTMG. Note that r = a corresponds 
to the boundary value of r, and that a Is also the tunnel 
radius. B0(K) needs to be loaded only If NBZDAT = 0. 

BILDP(K) Phase shift of voltage for cavity K; K = 1, LASTCV. 


v 


BCM(I) b 


BZDATA(I) 


DRDZ(I) dr 

ar 

FREQGH f 

I0BMA I 0 

INITCV 
ISAVE 


ISTATE 


I th beam radius In one-dimensional region, cm; 1=1, 

NUMB. See NUMB, NB1 , NB2 for related Inputs. If there Is 
no one-dimensional region (KAV13 = INITCV), only BCM(l) Is 
used and BCM(l) Is Initial beam radius. 

Experimentally obtained value of B z (r,z) at r = 0 and 
z = I th axial data point; 1=1, NBZDAT +1. 1=1 

corresponds to the beginning of a magnetic section and 
I = NBZDAT + 1 corresponds to the end of a magnetic 
section. The axial data points must be equally spaced. 
BZDATA(I) needs to be loaded only If NBZDAT f 0. 

Initial value of dr/d z for I th ring; j = i NRINGS. 

I = 1 corresponds to an Innermost ring. I = NRINGS corre- 
sponds to an outermost ring. 

Frequency, GHz 

Beam current at tube entrance, mA 
Number of first cavity considered In the presept case 
ISAVE = II + 12 + 13, where II, 12, and 13 ape defined as 
follows: 

II = 0: No action. 

11 = 1: Initial state of the present case will be stored 

In COMMON/STATE/for use In a future case. 

12 = 0: No action 

12 = 2: Initial state of the present case will be stored 

In a specified data set. 

13 = 0: No action. 

13 = 4: Final state of the present case will be stored 

In a specified data set. 

Used only when INITCV * 1 : 

ISTATE = 1: Initial state of present case Is equal to 

Initial state of previous case. 

ISTATE = 2: Initial state of present case Is equal to 

final state of previous case. 

ISTATE = 3: Initial state of present case Is equal to 

Initial state of the last case for which II (see ISAVE) 
was equal to 1 . 

ISTATE = 4: Initial state of present case Is equal to the 

state stored In a specified data set. 
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JSCF 

KAV13 

KIHP 

KLHAG 

KLOSS 

KPLOT 

KPPM 

KPRINT 

KREL 

KSHSIG 

KSPACE 

KSOLEN 

KVEL(I) 

KWRITM 

KWRITV 


LASTCV 

LASTMG 

LCIRCM(K) 

LGAPCH(K) 

LHAGCH(K) 

HSHAPE(K) 


Initially, space charge forces are evaluated once every 
OSCF steps. 

Number of the first cavity In the three-dimensional region 
KIHP = 0: Pierce Impedance Is entered as Input. 

KIHP = 1: Total Impedance Is entered as Input. (See ZIHP 

for related Input.) 

KLHAG * 0: Cavity lengths and magnetic sections coincide. 
KLHAG * 0: Cavity lengths and magnetic sections do not 
coincide. 


KLOSS = 0: ALPHLR(K) will be set equal to ALPHL(K) ; thus, 

ALPHLR(K) does not need to be loaded. 

KPLOT = 0: Plots are desired. 

KPLOT * 0: Plots are not desired. 

KPPH = 1: Single period PPH focusing Is to be used. 

KPPH = 2: Double period PPH focusing Is to be used. 

KPRINT = 0: Print Input data. 

KPRINT * 0: Do not print Input data. 

KREL = 0: Use relativistic equations of motion. 

KREL jt 0: Use nonrelatlvlstlc equations of motion. 

KSHSIG = 0: Print small-signal parameters. 

KSHSIG * 0: Do not print small-signal parameters. 

KSPACE = 0: Calculate space charge forces. 

KSPACE * 0: Set space charge forces equal to zero. 

KSOLEN = 0: Use PPH focusing. 

KSOLEN jt 0: Use solenoldal focusing; B z = BO, B r = 0. 

Number of cavity for the I*h printout of normalized axial 
velocities of rings; 1=1, NVEL 
KWRITH = 0: Do not print normalized masses of rings. 

KWRITH * 0: Print normalized masses of rings whenever 

normalized axial velocities of rings are printed. 

KWRITV = 0: Do not print normalized radial velocities of 

rings. 


KWRITV * 0: Print normalized radial velocities of rings, 

v r /u 0 , whenever normalized axial velocities of rings are 
printed. 

Number of last cavity considered In the present case 
( LASTCV < 70) 

Number of magnetic sections ( LASTHG < 70) 

Length of K th cavity, cm; K = 1 , LASTCV. 

Length of gap In K th cavity, cm; K = 1, LASTCV. 

Length of K th magnetic section, cm; K = 1, LASTHG. 
LHAGCH(K) needs to be loaded only If KLHAG £ 0. 

Used In defining electric field shape for K™ cavity, 
Inverse meters; K = 1, LASTCV; 


E z (a,z) at gap edge 
E z (a,z) at gap center 


= cosh (ml) 


where m = HSHAPE(K) and 1 = half-length In meters of 
K th gap. 
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MU( K) 


P 


Used In defining magnetic field shape for magnetic 

section. Inverse meters; K = 1, LASTMG; 


B z (a,z) at gap edge 
B z (a,z) at gap center = cos h(vg) 


NB1 

NB2 

NBWM 

NBZDAT 


NCAVSS 

NDISKS N d 

NHAX 

NPGRIO N r 

NPSC N sc ,r 

NRINGS R 

NUMB 

NVEL 

NXGRD1 N z 

NXGRD3 N z 

NXMAG 


where v = MU(K) and g = half-length In meters of K th 
magnetic gap. MU(K) needs to be loaded only If 
NBZDAT = 0. 

Used only when there Is a one-dimensional region 

(KAV13 > INITCV). Beam radius changes from first value to 
second value at end of NB1 th cavity. (NB1 needs to be 
loaded only If NUMB >1.) 

Used only when there Is a one-dimensional region 
(KAV13 > INITCV). Beam radius changes from second value 
to third value at end of NB2 th cavity. (NB2 needs to be 
loaded only If NUMB =3.) 

NBWM =0: No backward wave. 

NBWM > 0: Number of last cavity considered In calculating 

backward wave. 

NBZDAT = 0: Magnetic fields are determined by Input data 

B0 and MU. 

NBZDAT * 0: Magnetic fields are determined from experi- 

mental data. The number of experimental data values for 
B z (0,z) Is NBZDAT + 1. (NBZDAT < 100) 

Number of cavity for which small-signal parameters are 
calculated 

Number of disks (NDISKS < 24) 

For sums that are stored In tables, excluding space-charge 
force tables and magnetic field tables, NMAX Is the upper 
limit. (NMAX < 50) 

Number of grid points In r direction Is NPGRID +1. r = 0 
Is first grid point and r = a Is last grid point. 

Refers only to the grid for calculating forces on rings 
due to electric and magnetic fields. (NPGRID < 20) 

Number of grid points In r direction Is NPSC +1. r = 0 Is 
first grid point and r = a Is last grid point. Refers 
only to grid for calculating forces on rings due to 
space-charge fields. (NPSC < 20) 

Number of rings per disk (NRINGS < 4) 

Number of discrete values for the beam diameter In the one- 
dimensional region. 1 < NUMB < 3 

Number of cavities for which ring axial velocities are to be 
printed out. 

Number of nodes In z-dlrectlon per cavity. In one-dimensional 
region. (NXGRD1 < 64); NXGRD1 must be a multiple of 4. 

Number of nodes In z-dlrectlon per cavity. In three- 
dimensional region. (NXGRD3 < 64); NXGRD3 must be a 
multiple of 4. 

Number of grid points In z-dlrectlon per magnetic section. 
NXMAG must be a multiple of 4. Normally, NXMAG should be 
about the same as NXGRD3. NXMAG does not have to be 
loaded If KLMAG = 0. (NXMAG < 64) 
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PCOPA 

PINDBM 

TOLOV 


^ c /vj» a Normalized cathode flux: PCDPA = $ c /»|i a , where Is 

cathode flux at r = a. t|/ a = ira^B 0 , where a Is tunnel 
radius and B 0 = B0(1 ) . 

Input power, dBm 

Error criterion for determining whether an additional pass 
through a cavity Is required. If 


AV f.new ~ &V f.old 
u 

f .new 


> TOLOV 


TOLSC 


TOLTBL 


TWOGCM(K) 2g 
V0B V 0 

VJUMP(K) V Jk 
ZIMP(K) Z k 


then an additional pass Is made. 

Error criterion for determining whether the frequency of 
calculating space charge forces should be doubled. If 
for any ring the difference between the calculated space 
charge force and the extrapolation of the quadratic curve 
fit Is greater than the product of TOLSC and the maximum 
rf electric field force In the cavity, the frequency of 
calculating space charge forces Is doubled. 

Error criterion used In electric and magnetic field tables 
for determining the number of terms to Include In Infinite 
summations. Suggested value: TOLTBL = 0.02. 

Length of gap In magnetic section, cm; K = 1, LASTMG 

Beam voltage, volts 

The dc voltage jump for cavity, volts; K = 1, LASTCV 

If KIMP = 0, ZIMP(K) Is the Pierce Impedence In ohms for 
K th cavity; K = 1, LASTCV 

If KIMP = 1, ZIMP(K) Is the total Impedance In ohms for 
K th cavity; K = 1, LASTCV 


DESCRIPTION OF OUTPUT DATA 

The program output consists of both printed output and plots. The plots 
are produced only If the Input parameter KPLOT Is zero. The printed output 
Includes three parts: (1) printing of small-signal parameters (only If the 

Input parameter KSMSIG Is zero); (2) a cavlty-by-cavlty printing of selected 
data; (3) printing of data on ring dynamics at selected cavities. The program 
can be easily modified to calculate and print other data If desired. 

First, the Input data are printed In NAMELIST format. (The printing of 
the Input data can be suppressed by letting KPRINT be nonzero.) Then the 
small-signal parameters are printed If desired. The small-signal parameters 
are as follows: 

Name Symbol 

U0 u 0 

BEB 0 e b 

BIB Bib 

KP K p 

ZC Z c 

C C 

B b 

D d 


Description 

Initial beam velocity, m/sec 

Product of beam propagation constant and beam radius 

Product of wave propagation constant and beam radius 

Pierce Impedance, ohms 

Total Impedance, ohms 

Pierce's C 

Pierce's b 

Pierce loss parameter d 
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DGAIN 

QQ 

A1PA2 


Asymptotic gain per cavity, dB/cavlty 

Q Pierce's space charge parameter, QC/C 

A] + Launching loss, dB 


The following data are printed for each cavity: 


CAV 

VMAG 

I SHAG 

ISPHA 

LSGAIN 

POUT 

AVERHO 

RMSANG 


Cavity number 

Magnitude of gap voltage, volts 

Magnitude of normalized Induced current 

Phase of normalized Induced current divided by ir 

Large-signal gain, dB 

Output power divided by I 0 V 0 

-» NK 

Average of the ring radii P = — ^ p^ (NR = NDISKS*NRPD) 

RMS value of the angle that the velocity vectors of the rings 
make with the z-axls; In degrees 


RMSVEL 


PKE 

INTRC 

PRF 

PBW 

PLC 

PBAL 


SC 


RMSANG = 



180 

ir 


RMS value of the normalized radial velocities of the rings, 
where the normalization Is with respect to u n , the 
Initial axial velocity 



Change In beam kinetic power referred to Initial beam kinetic 
power, divided by I 0 V 0 

Power loss due to beam Interception, divided by I Q V n 
Power In forward wave, divided by I 0 V 0 
Power In backward wave, divided by I 0 V 0 
Cumulative power loss (except for loss due to beam Inter- 
ception), divided by I 0 V 0 
Power balance equal to 


1 + PKE - PJUMP + PRF + PBW + PLC + INTRC 


where PJUMP Is the power due to all voltage jumps up to the 
current cavity 

Integer variable KSCF: space charge forces are evaluated 

every KSCF Integration steps. 
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At selected cavities, the following data are printed: 

"Normalized axial velocities" - Axial velocities normalized with respect to 

u 0 , the Initial axial velocity 

"Normalized radial velocities" - Radial velocities normalized with respect to 

u 0 , the Initial axial velocity; printed 
only If KWRITV * 0 

"Normalized ring masses" - Ring masses normalized with respect to Initial 

ring mass; printed only If KWRITM * 0 

"Ring radius divided by a" - Ring radii normalized with respect to a 


DESCRIPTION OF GLOBAL VARIABLES 


Global variables are variables that are used by more than one subroutine 
and hence are put Into COMMON blocks. In the following descriptions, we ex- 
clude Input data variables since these have already been described. 


Name 

Symbol 

DELV(K) 

AVfk 

DELVR(K) 

AV bk 

OVTEMP(K) 

AVfk 

DVRTEM(K) 

AVbk 

VGAP(K) 

Vk 

VGAPF(K) 

Vfk 

VGAPR(K) 

Vbk 

A 

ALFI(J.I) 
ALF2( 0,1) 
ALF10(J) 

a 

ALF20(J) 

ALPHLB(K) 

(“Dbk 

B( I) 

b 


Description 

Forward Induced voltage for k*h cavity, 
volts; K = 1, LASTCV 

Backward Induced voltage for cavity, 

volts; K = 1, LASTCV 

Approximate forward Induced voltage for k^^ 1 
cavity, obtained by assuming rings have 
constant axial and radial velocities 
throughout cavity, volts; K = 1, LASTCV 
Backward Induced voltage of previous case, 
for kth cavity, volts; K = 1, LASTCV 
Sum of forward voltage and backward voltage 
for k th cavity, volts; K = 1, LASTCV 
Forward voltage for k** 1 cavity, volts; 

K = 1, LASTCV 

Backward voltage for k** 1 cavity, volts; 

K = 1, LASTCV 
Tunnel radius, meters 
Ratio of Inner radius to mean radius of 
(J,I) th ring, 1=1, NDISKS; 0=1, NRINGS 
Ratio of outer radius to mean radius of 
(J,I) th ring, 1=1, NDISKS; J = 1, NRINGS 
Initial ratio of Inner radius to mean radius 
of a ring, where J = 1 Is an Innermost 
ring and J = NRINGS Is an outermost ring 
Initial ratio of outer radius to mean radius 
of a ring, where J = 1 Is an Innermost 
ring and J = NRINGS Is an outermost ring 
Voltage attenuation of backward wave In going 
from cavity K + 1 to cavity K, dB per 
cavity; K * 1, LASTCV 
I™ beam radius In one-dimensional region, 
meters; 1=1, NUMB 
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i 

j 


BSUBR(IR.IZ) 

BSUBZ(IR.IZ) 

| 

i 

C 

C1(I) 

C2(I) 

CAVFLG 


CNORM 

CONST1 

CONST2 

C0NST4 

COSY(K) 

COSYB(K) 


CTHETA 

DELP 


OELPSC(I) 


OELX(I) 

DELXSC(I) 


DPR 

ECHARG 

EMASS 

EPS0 

ERRWT( I , J) 


FREQ 
1 0B 


B r( r IR* z Iz) Radial component of focusing magnetic field 
at grid point webers per square 

meter; IR = 1 , NPGRID +1, IZ = 1, NXMAG. 
Calculated In TBL3H and MAGCF. 

B z( r IR* z Iz) Ax 1 a "l component of focusing magnetic field 
at grid point (rjp.Zj^), webers per square 
meter; IR = 1, NPGRID +1, IZ = 1, NXMAG. 
Calculated In TBL3M and MAGCF. 

c Speed of light, meters per second 

Constant associated with I th beam diameter; 
1=1, NUMB. Calculated In INDAT. 

Constant associated with I th beam diameter; 
1=1, NUMB. Calculated In INDAT. 

Constant calculated In INDAT. 

Flag used by subroutine KSCKC to determine 
which of two adjacent cavities Is to be 
selected when positioned at the boundary 
between the two cavities 

c/ua Normalized speed of light 

Constant calculated In INDAT 

No longer used 

Constant calculated In INDAT 

Array of cosine values used In TBL3E; 

K = 1, 2 * NXGRID + 1. Calculated In INDAT. 

Array of cosine values used In TBL3M; 

K = 1, 4 * KPPM * NXMAG + 1. Calculated In 
INDAT. 

cos 0 Cosine of normalized time 

Normalized distance between two adjacent grid 
points In radial direction. Calculated In 
INDAT. 

I th normalized distance between two adjacent 
grid points In radial direction In space 
charge tables; 1=1,2. Calculated In 
SCAT3. 

Integration step size In I th cavity string; 
1=1, NSCAV. Calculated In INDAT. 

I th normalized distance between two adjacent 
grid points In axial direction In space 
charge tables; 1 = 1 , 2 . Calculated In 
SCAT3. 

Conversion factor from degrees to radians, 

DPR = 1 80/ir 

e Charge of an electron, coulomb; ECHARG Is 

positive 

Mass of an electron, kg 

c o Permittivity of free space, farad/meter 

Error weights used In determining the number 
of terms to Include In the summations 
calculated In TBL3E and TBL3M; 1=1, 5; 

J = 1, NPGRID + 1. Calculated In INDAT. 

f Frequency, Hz 

Io Beam current at tube entrance, amps 
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I0BAVE 


Beam current at tube entrance, reduced by the 
percentage of the beam that has been Inter- 
cepted, at the present axial position, amps 
Calculated In ACC. 

INTLS 


Power loss due to the cumulative beam Inter- 
ception that has taken place up to the pre- 
sent axial position, watts 

ISCAV(K) 


Number of cavity string In which K th cavity 
Is located; K = 1, LASTCV 

ITER 


Present number of Iteration of beam rings 
traversing cavity 

KBEGIN 


Equals zero at beginning of case and Is set 
equal to 1 Immediately thereafter. 

KCAV 


Number of current cavity 

KOUNTP 


Counter for data points to be plotted 

KPOZC(K) 

Kp/Z 

Ratio of Pierce Impendence to total Impedance 
for K th cavity; K = 1, LASTCV 

KSCAV 


Number of present cavity string. (A cavity 
string Is a string of consecutive cavities 
all having the same geometrical and elec- 
trical properties.) 

KSCF 


Space charge forces are calculated every 
KSCF steps. 

KSCFT 


Used to store a previous value of KSCF. 

KSTOP 


Equals zero until end of last cavity, when 
It Is set equal to 1 . 

KX 


Number of current node, counted from first 
node of first cavity 

KX0 


Initial node of present case, counted from 
first node of first cavity 

KXPRI 


Number of current node, counted from first 
node of current cavity 

KXSC 


Number of node, counted from first node of 
first cavity, at which the most recent cal- 
culation of space charge forces was done 

LAMB 

X 

Beam wavelength, meters 

Length of K* h cavity, meters; K = 1, LASTCV 

LCIR(K) 

Lk 

LDISK 

*d 

Axial length of a ring, meters 

LDISKN 

H d /a 

Normalized axial length of a ring 

LMAG(K) 

Length of K** 1 magnetic section, meters; 
K = 1, LASTMG 

LSMALL(K) 

lk 

One-half of gap length In K** 1 cavity, 
meters; K = 1, LASTCV 

MDISK 

("d 

Mass of a disk at tube entrance, kg 

MLSHML 

mi 

mi Is the product of MU and LSMALL for the 

MRING0 

slnh mi 

current cavity 


Mass of a ring at tube entrance, kg 
Mass of (J,I) th ring, kg; 1=1, NOISKS; 
J = 1, NRINGS 

MRING(J.I) 


MUA(I) 


I*h Z ero of Bessel function J 0 ; I = 1 100 

MU0 

>*0 

Permeability of free space, henry/meter 

NB 


Number of current beam-diameter region. 

(A beam-diameter region Is a region 
throughout which the beam diameter Is con- 
stant.) Applies only to one-dimensional 
region. 
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NCAVPS(I) 

NCASE 

NCUMCS(I) 

NPSC1 

NR 

NRPD 

NSCAV 

NXGRID 

NXISC 

POB 

PBW(K) 

PHIDN(J.I) 


PI 

PKE(K) 

PLC 

POBAL 

PRF(K) 

PRICL 


PSI(IR.IZ) 


PSIA 

PSC1(J,I) 
PSC2( J , I) 
PSC3( J , I) 


PSCA( J,I) 
PSCB(J.I) 
PSCC(J.I) 

QDISK 

QRINGO 

QRING(J.I) 

RELC(J.I) 


P o 




If 


»H r IR,Zlz) 


wa 2 B 0 


<ld 



Number of cavities In I*h cavity string; 
1=1, NSCAV 

Number of current case 

Number of last cavity In I th cavity string; 
1=1, NSCAV 

Number of grid points In the first of two 
radial grids In space-charge tables 

Total number of rings In a beam wavelength 
( NDISKS * NRPD) 

Number of rings per disk (NRPD = 1 In one- 
dimensional region.) 

Total number of cavity strings 

Number of nodes In z-dlrectlon per cavity 

Number of grid points In the axial grid In 
space-charge tables 

I 0 /u 0 , charge density per unit length at 
tube entrance, coulomb/m 

Power In backward wave In K** 1 cavity, 
watts; K = 1, LASTCV. Calculated In ENBAL. 

Normalized angular velocity of (J-I) th ring; 
1=1, NDISKS; J = 1, NRINGS. Calculated In 
FACC3. 

3.141593 

Kinetic power at end of K^ cavity minus 
kinetic power at beginning of tube, watts; 

K = 1, LASTCV. Calculated In ENBAL. 

Cumulative power loss at end of K th cavity, 
watts; K = 1, LASTCV. Calculated In ENBAL. 

Power balance. Calculated In ENBAL. 

Power In forward wave In K th cavity, watts; 

K = 1, LASTCV. Calculated In ENBAL. 

Power In backward wave In present cavity, 
watts. Same as PBW(K). Calculated In 
ENBAL. 


Magnetic flux at grid point (rjp,ZT 7 ), webers; 
IR = 1, NPGRID +1, IZ = 1, NXMAG. 

Calculated In TBL3M and MAGCF. 

Normalizing factor. Calculated In INDAT. 
Radial space charge force on (J,I)th ring 
at each of the latest 3 selected values of 
Independent variable \ used for quadratic 
curve-fitting; J = 1, NRINGS; 1=1, NDISKS. 
Calculated In SCFT3. 


Coefficients of the quadratic polynomial used 
for curve-fitting the radial space charge 
force on (J,I) th ring; J = 1, NRINGS; 
1=1, NDISKS. Calculated In SCFT3. 

Charge of a disk at tube entrance, coulombs 
Charge of a ring at tube entrance, coulombs 
Charge of (J-I)™ ring, coulomb; 

1=1, NDISKS; J = 1, NRINGS 


Relativistic correction factor for (J-I)th 
ring; 1=1, NDISKS; J = 1, NRINGS. ' 
Calculated In FACC3. 
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RHO0( J) 


Initial value of normalized mean radius of a 
ring, where J = 1 Is an innermost ring and 
J = NRINGS Is an outermost ring; 

J = 1 , NRINGS. Calculated In INDAT. 

RHOMAX 


Normalized outer radius of the outermost ring 
at the current value of Calculated 

In ACC. 

RHO( J,I) 


Normalized radius of (J-I)th ring; 

1=1, NDISKS; J = 1, NRINGS. Calculated 
In XINT. 

RHOP( J,I) 


First derivative of p^ with respect to ?; 

RHOPP(J.I) 

dE 

1=1, NDISKS, J = 1 , NRINGS. Calculated 
In XINT. 


Second derivative of p^ with respect to 5 ; 



1=1, NDISKS, J = 1, NRINGS. Calculated 
In ACC. 

SBET0A(I) 


Value of B 0 a for cavities In I ** 1 cavity 
string; 1=1, NSCAV 

SCA1(I) 


ith entry In table of space charge forces 

SCA2( I) 


for first, second, and third beam diameters; 

SCA3 ( I) 


1=1, NXISC + 1. Calculated In SCAT1 . 
Applies to one-dimensional region only. 

SCHOLD 


No longer used. 

SIMPS(I) 

*1 

ith coefficient for Simpson's rule Integra- 
tion; 1=1, NXGRID. Calculated In INDAT. 

SINY(K) 


Array of sine values used In TBL3E; 

K = 1, 2 * NXGRID + 1. Calculated In INDAT. 

SINYB(K) 


Array of sine values used In TBL3M; 

K = 1, 4 * KPPM*NXMAG + 1. Calculated In 
INDAT. 

SLCIR(I) 


Length of cavities In 1^ cavity string, 
meters; 1=1, NSCAV 

SLSMAL(I) 


One-half of gap length for cavities In I** 1 
cavity string, meters; 1=1, NSCAV 

SMSHAP(I) 


Electric field shaping factor for cavities In 
ith cavity string, meters; 1=1, NSCAV 

SR1MB2 

/ 7T 

.2 Relativistic correction factor. Calculated 


Vi-G* 

J In INDAT. 

SSGAIN(K) 


/ 

Small-signal gain at K th cavity, based on 
small-signal calculations using geometrical 
and electrical properties of cavity number 
NCAVSS; K = 1, LASTCV. Calculated In 
SMSIG. 

STHETA 

sin 6 

Sine of normalized time 

THT( J,I) 


Value of normalized time of arrival of 
(J-I) th ring at current axial position 
H 1=1, NDISKS; J = 1, NRINGS. 
Calculated In XINT. 

THTP(J,I) 

de ii 

d* 

First derivative of e., with respect to 5 ; 
1=1, NDISKS; J = 1, NRINGS. Calculated 
In XINT. 
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THTPP(J.I) 

“St 


d ? 2 

THTPLT(J.I) 


TPXR( IR, IZ,L) 
TPXI( IR, IZ,L) 


TXR( IZ, L) 
TXI(IZ.L) 


TWOG(K) 

2g 

TWOPI 

2ir 

U0B 

Uq 

U0NORM 

V0BSUM(K) 

u 0 /“a 

VGAP1M 

VMAX 


W 

G) 

X 

S 

X0 

So 


XICAV(K) 

XIMAG(K) 

XISCAV(I) 

XIHXSC 


XISC1 

XISC2 

XISC3 

XPRI 


XSC1(J,I) 
XSC2( J,I) 


Second derivative of e u with respect to 5; 
1=1, NDISKS: J = 1 , WRINGS. Calculated 
In ACC. 

Value of e for (J,I) th ring used In plots; 

1=1, NDISKS; J = 1, NRINGS 
Real (TPXR) and Imaginary (TPXI) components 
of the L™ table for the data point 
( r IR* z IZ ) • IR = 1. NPGRID + 1; 

IZ = 1, NXGRID: L = 1 , 5 for TPXR: L = 1, 3 
for TPXI. Calculated In TBL3E. 

Real (TXR) and Imaginary (TXI) components of 
the L tf1 table for the data point ZT7; 

IZ = 1, NXGRID; L = 1, 3 for TXR: 1=1,2 
for TXI. Calculated In TABL1 . Used for 
one-dimensional region only. 

Gap length of magnetic section; meters; 

K = 1, LASTMG 
2ir 

Initial axial beam velocity, m/sec 
Normalized Initial axial beam velocity 
Beam voltage plus sum of all voltage jumps up 
to K th cavity, volts; K = 1, LASTCV 
Voltage across first gap, volts 
The maximum |VGAP(K)|, where the maximum Is 
taken over all values of K up to the pre- 
sent cavity; volts. Calculated In BTWNC. 
Radian frequency, u = 2wf; rad/sec 
Current value of normalized position along 
tube axis, measured from beginning of first 
cavity 

Initial value of normalized position along 
tube axis for the present case, measured 
from beginning of first cavity 
Normalized axial position at end of K th 
cavity; K = 1, LASTCV 
Normalized axial position at end of 
magnetic section; K = 1, LASTMG 
Normalized axial position at end of I th 
cavity string; 1=1, NSCAV 
Maximum normalized axial separation distance 
In space charge force tables. (If two 
rings are separated by a distance greater 
than XIMXSC, the space charge force on one 
ring due to the other Is set to zero.) 
Calculated In INDAT. 

The three latest values of 5 at which space 
charge forces were calculated 

Present value of normalized position along 
tube axis, measured from beginning of cur- 
rent cavity. Calculated In ACC. 

Axial space charge force on (J,I) th ring at 
each of the latest 3 selected values of In- 
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XSC3( J,I) 


XSCA( J,I) 

XSCB( J , I) 

XSCC( J , I) 

ZC(K) 

IOTABL(K) I 0 (x k ) 

IITABL(K) I^Xk) 

SCAP( IP1 , IP2, IX, II) 

SCAX( IP1 , IP2, IX, II) 


dependent variable 5 used for quadratic 
curve-fitting; J = 1, NRINGS, 1=1, NDISKS. 
Calculated In SCFT3. 

Coefficients of the quadratic polynomial used 
for curve-fitting the axial space charge 
force on (J-I) th ring; J = 1, NRINGS; 

1=1, NDISKS. Calculated in SCFT3. 

Total Impedance for cavity, ohms; 

K = 1, LASTCV 

Tables of Bessel functions I 0 and I-j, 
starting with X] = 0 and with an Incre- 
ment of 0.03 In x; K = 1,5667 

Tables of radial (SCAP) and axial (SCAX) 
space charge forces for three-dimensional 
region; IP1 =1, NPSC + 1; IP2 = 1, NPSC + 1; 
IX = 1, NXISC +1; II = 1, NRINGS. Calcu- 
lated In SCAT3. 


DESCRIPTION OF SUBROUTINES 

In this section, the COMMON blocks, the main program, and each subroutine 
are described. Flow charts are Included for the more complex subroutines. 


COMMON Blocks 

The program has five COMMON blocks. COMMON/INDATA/ contains the Input 
data. COMMON/MREAL/ contains real and Integer global variables. COMMON/MCOMP/ 
contains complex global variables. COMMON/PLOT/ contains arrays used for 
plotting Applegate diagrams. COMMON/STATE/ contains arrays used for storing 
Initial states and final states. An Initial state Is a collection of vari- 
ables whose values are those at the beginning of the first cavity In the tube 
section under consideration. A final state Is the same collection of vari- 
ables but whose values are those at the end of the last cavity In the tube 
section. The variables In this collection are those whose values must be known 
In order to start the simulation of a tube section. It Is necessary to store 
these Initial and final states because they may be used In later cases. For 
example. In simulating the second tube section, we use the final state of the 
first tube section as a startup. 


Main Program 

The main program reads the first 100 zeros of the Bessel function Jn. 
Subroutines STDAT, INDAT, and SMSIG are then called. When the last Integration 
step has been done, the final state Is stored In COMMON/STATE/. Depending on 
the value of ISTATE, the final state may also be stored In a data set for 
possible use In a future case. The main program calls XINT to Integrate one 
step and calls BTWNS to check on what action should be taken between steps. 
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Subroutine ACC 


Subroutine ACC Is called by subroutine XINT. Subroutine ACC first checks 
to see If a new space-charge calculation Is required. A new calculation Is 
required If the last calculation was KSCF steps ago or If the current value of 
% Is at a boundary between two beam-diameter regions. If a new calculation 
Is required, then either subroutine SCFT1 or SCFT3 Is called. Next, subroutine 
ACC calculates the normalized axial and radial accelerations of each ring. In 
the one-dimensional mode this calculation Is done In the same way as In the 
one-dimensional computer program (ref. 4). In the three-dimensional mode, the 
accelerations result from the following forces: space-charge forces, rf 

electric-field forces, voltage-jump forces, magnetic forces, centripetal forces, 
and forces due to relativistic effects. The space-charge forces are calculated 
by quadratic curve fits, and the other forces are obtained by a call to sub- 
routine FACC3. Subroutine ACC also calculates THTPP(J,I) and RH0PP(J,I), the 
second derivatives of the functions OijU) and PijU), respectively. 


Subroutine APPLE 

Subroutine APPLE Is called by the main program at the end of the simula- 
tion. Subroutine APPLE manipulates plot data that has been stored In the 
course of the simulation and calls the plotting device for plotting r - z 
plots and, optionally, t - z plots. 


Subroutine BTOS 

Subroutine BTOS Is called by subroutine INDAT and by the main program. 
When called by INDAT, subroutine BTOS stores In COMMON/STATE/ the collection 
of variables making up the Initial state. When called by the main program, 
subroutine BTOS stores In COMMON/STATE/ the collection of variables making ud 
the final state. 


Subroutine BTWNC 

Subroutine BTWNC Is called by subroutine BTWNS when the end of a cavity Is 
reached. Subroutine BTWNC determines whether another pass through the cavity 
Is required. A second pass Is not required If either of two conditions Is met: 
(1) a second pass has already been made, or (2) the Induced forward voltage Is 
sufficiently close to the approximation of the Induced forward voltage obtained 
by assuming constant ring velocities throughout the cavity. 

If a second pass Is not required, the subroutine sets up for entrance Into 
the new cavity. If needed, new tables are calculated by calling the appro- 
priate subroutines. Since a second pass through the new cavity may be needed, 
It Is necessary to store the state at the beginning of the new cavity. Sub- 
routine FSAPP Is called to calculate DVTEMP, the approximation to the Induced 
forward voltage for the new cavity. The new forward voltage Is then calcu- 
lated by phase shifting and attenuating the old forward voltage and adding 
DVTEMP. If the old cavity was the last cavity In the one-dimensional region 
and the new cavity Is the first cavity In the three-dimensional region, sub- 
routine TRN13 Is called. 
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Subroutine BTWNS 


Subroutine BTWNS Is called by MAIN. Subroutine BTWNS determines what 
action should be taken between Integration steps. First, subroutine INTRC Is 
called to check on beam Interception. Next, the subroutine determines whether 
the end of a cavity has been reached. If so, subroutine BTWNC Is called. If 
the end of a magnetic section has been reached, subroutine TBL3M Is called. 
Finally, the subroutine determines whether subroutine PLOTS should be called. 
Subroutine PLOTS Is called only In the three-dimensional region. 


Subroutine CAVP 

Subroutine CAVP Is called by subroutine INDAT. Subroutine CAVP calculates 
quantities related to cavities and cavity strings. A cavity string Is a string 
of consecutive cavities all having the same electrical and geometrical proper- 
ties. The number of cavities cannot exceed 70, and the number of cavity 
strings cannot exceed 35. 


Subroutine CKPDZ 

Subroutine CKPDZ Is called by subroutine INDAT. Subroutine CKPDZ calcu- 
lates KPDZC(K) , the ratio of Pierce Impedance to total Impedance, for the K th 
cavity. 


Subroutine CR00T 

Subroutine CROOT Is called by subroutine SMSIG. Subroutine CROOT calcu- 
lates the roots of the quartlc polynomial Involved In calculating small-signal 
parameters. The subroutine first solves for the roots of the cubic polynomial 
using the Newton-Raphson method. Using these roots as the first Iterates, the 
subroutine then uses the Newton-Raphson method to solve for the roots of the 
quartlc polynomial . 


Subroutine DVSUM 

Subroutine DVSUM Is called by subroutine XINT. In the three-dimensional 
region, DVSUM Is also called by subroutine FSAPP. Subroutine DVSUM calculates 
one term In the summation for the Induced forward voltage, DELV, and one term 
In the summation for the Induced backward voltage, DELVR. These terms are 
added to the running sums for DELV and DELVR. The summations are those given 
In equations (16) and (17). 


Subroutine ENBAL 

Subroutine ENBAL Is called by subroutine BTWNC. Subroutine ENBAL calcu- 
lates the power balance at the end of the cavity In question. There are six 
terms In the power balance: PKE, PJUMP, PRF, PBW, PLC, and INTLS. PKE Is the 

kinetic power minus the Initial kinetic power. PJUMP Is the power due to all 
voltage jumps up to the cavity In question. PRF Is the power In the forward 
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wave minus the power In the forward wave In the first cavity. PBW Is the power 
In the backward wave In the first cavity minus the power In the backward wave. 
PLC Is the sum of power losses up to the present cavity. INTLS Is the power 
loss due to beam Interception. Power balance Is equal to one plus the sum of 
these six terms divided by the Initial beam power. 


Subroutine FACC1 

Subroutine FACC1 Is called by subroutine ACC when the simulation Is In 
the one-dimensional mode. Subroutine FACC1 calculates EXACC, the normalized 
acceleration of a ring due to rf electric-field force and voltage-jump force. 


Subroutine FACC3 

Subroutine FACC3 Is called by subroutine ACC when the simulation Is In 
the three-dimensional mode. Subroutine FACC3 calculates EXACC and EPACC, the 
normalized axial and radial accelerations of a ring due to rf electric-field 
force and voltage-jump force. The subroutine then calculates MXACC and MPACC, 
the normalized axial and radial accelerations due to magnetic focusing force. 
The normalized centripetal acceleration, CENTRP, Is then calculated. Finally, 
the subroutine calculates RXACC and RPACC, the normalized axial and radial 
accelerations due to relativistic effects. 


Subroutine FIELD 

Subroutine FIELD Is called by subroutine FACC3 when the simulation Is In 
the three-dimensional mode. Subroutine FIELD calculates EXRO, EXIO, EPRO, 
EPIO, BXO, BPO, PSIO, EXDCV, and EPDCV. EXRO and EXIO are the real and 
Imaginary parts of the shaping factor for the axial electric field. EPRO and 
EPIO are the real and Imaginary parts of the shaping factor for the radial 
electric field. BXO, BPO, and PSIO are the shaping factors for the axial 
magnetic field, radial magnetic field, and magnetic flux field, respectively. 
EXDCV and EPDCV are the shaping factors for the axial and radial voltage-jump 
fields. 


Subroutine FSAPP 

Subroutine FSAPP Is called by subroutine BTWNC . Subroutine FSAPP calcu- 
lates an approximation, DVTEMP, to the Induced forward voltage by assuming 
the rings have* constant axial and radial velocities In the cavity under 
consideration. 


Subroutine INDAT 

Subroutine INDAT Is called by the main program. The subroutine first 
reads Input data from a prestored data set. It then reads data entered at the 
terminal In NAMELIST format. Only Input data that differ from those In the 
prestored data set need to be entered. The Input data are stored In a data 
set for possible future use. The subroutine then calculates constants and 
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Initializes variables. The Initial state Is obtained In one of three ways: 

(1) If the Input parameter ISTATE Is 4, the Initial state Is read from a data 
set; (2) If ISTATE Is not 4 and INITCV (Initial cavity) Is 1, subroutine INIT 
Is called and the Initial state Is calculated from Input data; (3) If ISTATE 
Is not 4 and INITCV Is not 1, subroutine STOB Is called and the Initial state 
Is obtained from COMMON/STATE/. There Is storage allotment for three states 
In COMMON/STATE/. Which of the three Is chosen Is determined by ISTATE (see 
definition of ISTATE In the section on description of Input data). 

After the Initial state Is obtained. It Is stored In COMMON/STATE/ In the 
first of the three storage allotments. This Is done because the Initial state 
of the present case may be used as the Initial state In the following case, 
and thus It must be stored. The Initial state Is also stored In the third of 
the three storage allotments If the Input parameter ISAVE Is an odd Integer. 
This Is done because the Initial state of the present case may also be required 
as the Initial state In some future case after the following case, and thus It 
must be stored. 

The Input parameter ISAVE Is the sum of II, 12, and 13 (see definition of 
ISAVE In the section on description of Input data). If 12 Is 2, the Initial 
state of the present case Is stored In a specific data set for possible use In 
a future case. 

If INITCV Is 1, subroutine INDAT calculates the forward voltage for the 
first cavity. If INITCV Is not 1, the forward voltage for the cavity previous 
to cavity INITCV Is required for startup. This voltage will be known either 
from the previous case or by reading It from a data set. The backward voltages 
are calculated next. Finally, the Induced backward voltages OELVR are stored 
In the array DVTREM. This Is done because In the equations of motion we must 
use the Induced backward voltages of the previous case. We cannot use the 
OELVR array because this array changes during the simulation as the new Induced 
backward voltages are calculated. We therefore use the DVRTEN array In the 
equations of motion. 


Subroutine INIT 

Subroutine INIT Is called by subroutine INDAT and only when INITCV (Ini- 
tial cavity) Is 1. Subroutine INIT calculates Initial values for the collec- 
tion of variables making up the Initial state. 


Subroutine INTRC 

At the end of each Integration step, subroutine BTWNS calls subroutine 
INTRC to check for beam Interception. The check Is done for each ring. One 
of four conditions may occur: (1) both the Inner and the outer radius of the 

ring are less than a, (2) the outer radius Is greater than a and the Inner 
radius Is less than a, (3) both the Inner and the outer radius are greater 
than a, (4) the ring has already been completely Intercepted. 

If condition (1) occurs, no action Is taken. If condition (2) occurs, 
the outer radius Is redefined to be a, and the portion of the ring located 
beyond r = a Is assumed to be Intercepted. The remaining partial ring, since 
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Its outer radius has been redefined, has a redefined value of centroid radius 
given by 


r = 


The mass and charge of the ring are reduced by the appropriate amounts. The 
power lost due to the Interception Is calculated and added to the running 
total. If condition (3) occurs, the ring Is assumed to be completely Inter- 
cepted. Instead of setting the mass and charge of the ring to zero. It Is 
more convenient to set these quantities to negligible amounts and keep the 
ring In the simulation. If condition (4) occurs, the outer radius of the 
Intercepted ring Is set to a, and the Inner radius Is set to a quantity 
slightly less than a. 



Subroutine MAGCF 

Subroutine MAGCF Is called by subroutine TBL3M. In subroutine TBL3M, the 
arrays BSUB2B, BSUBRB, and PSIB are calculated prior to calling subroutine 
MAGCF. These arrays correspond to the magnetic quantities B z , B r , and 
evaluated at the grid points of the magentlc section In question. The axial 
grid points are equally spaced over the length of the magnetic section. If 
each magnetic section does not coincide with a cavity, these axial grid points 
will not coincide with the Integration nodes. The program requires that the 
axial grid points for the B z , B r , and tables coincide with the Integra- 
tion nodes. Subroutine MAGCF calculates these tables by performing linear 
Interpolations on the BSUBZB, BSUBRB, and PSIB arrays. The results are stored 
In the BSUBZ, BSUBR, and PSI arrays. 


Subroutine NBEAM 

Subroutine NBEAM Is called by subroutines BTWNC and SCFT1 . When the 
simulation Is In the one-dimensional mode, subroutine NBEAM calculates the 
number of the beam-diameter region In which XXX Is located. 


Subroutine OUTPT 

Subroutine OUTPT Is called by subroutine BTWNC. Subroutine OUTPT calcu- 
lates and prints output data. Since the output data have already been de- 
scribed In another section, we will not repeat the description here. 


Subroutine PLOTS 

Subroutine PLOTS Is called by subroutine BTWNS at each grid point z for 
which plotting Information Is to be stored. Two types of plots are used: r 

as a function of z and t as a function of z, for each ring to be plotted. 
The latter type of plot, which Is similar to the Applegate diagram. Is 
optional. The values of z, r, and t for the rings In question are stored 
In arrays for later usage by the plotting device. 
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Subroutine SACC1 


Subroutine SACC1 Is called by subroutine SCFT1 (one-dimensional mode 
only). Subroutine SACC1 calculates the space-charge acceleration of one ring 
due to another by linear Interpolation on space-charge force tables. Although 
we use the term space-charge force tables, these tables actually evaluate 
accelerations rather than forces. If the two rings are In two different beam- 
diameter regions, two linear Interpolations are done and the average Is taken. 


Subroutine SACC3 

Subroutine SACC3 Is called by subroutine SCFT3 (three-dimensional mode 
only). Subroutine SACC3 calculates the axial and radial space-charge accelera- 
tion on one ring due to another by a three-dimensional linear Interpolation on 
space-charge force tables. Although we use the term space-charge force tables, 
these tables actually evaluate accelerations rather than forces. 


Subroutine SCAT1 

Subroutine SCAT1 Is called by subroutine INDAT or subroutine INIT. For 
each beam diameter IB and for the I** 1 separation distance X, subroutine 
SCAT1 calculates the space-charge acceleration SCATHP(I.IB) of one disk due to 
another for the one-dimensional mode. Using EQUIVALENCE statements, we put 
SCA1, SCA2, and SCA3 Into COMMON/MREAl/ Instead of SCATMP. The arrays SCA1 , 
SCA2, and SCA3 make up the space-charge force tables for the one-dimensional 
mode. In this subroutine, X Is used as the separation distance between disks. 
Since X Is also used as a global variable (current value of 5), It Is 
necessary to save the global variable In XSAVE. At the end of the subroutine, 
the value of the global variable Is returned to X. In addition, subroutine 
SCAT1 calculates the ratio of the largest space-charge acceleration to the 
smallest. The smallest space-charge acceleration occurs when the separation 
distance Is largest. Subroutine SCAT1 also calculates the ratio of the beam 
wavelength to the maximum separation distance. 


Subroutine SCAT3 

Subroutine SCAT3 Is called by subroutine INOAT or subroutine INIT. For 
the three-dimensional mode, subroutine SCAT3 calculates and stores In a table 
the axial and radial space-charge accelerations of one reference ring due to 
one source ring. The table Is three-dimensional since these accelerations are 
functions of three variables: (1) r-| , the centroid radius of the source ring; 

(2) P2, the centroid radius of the reference ring; (3) z, the axial separation 
between the rings. This subroutine first calculates the number of z grid 
points and the spacing between z grid points. There are two different grid 
spaclngs for z (DELXSC(I), I =1, 2). Similarly the subroutine calculates 
two different grid spaclngs for the r variables (DELPSC(I), I =1, 2). The 
subroutine proceeds to calculate the space-charge accelerations for every grid 
point. In addition, subroutine SCAT3 calculates the ratio of the largest 
space-charge acceleration to the smallest. The smallest space-charge accelera- 
tion occurs when the axial separation distance Is largest. Subroutine SCAT3 
also calculates the ratio of the beam wavelength to the maximum separation 
distance. 


47 


Subroutine SCFT1 


Subroutine SCFT1 Is called by subroutine ACC whenever a new calculation of 
space-charge accelerations In the one-dimensional mode Is required. Subroutine 
SCFT1 first calculates what the space-charge accelerations would be If the pre- 
sent curve fits are used. Then for each disk the subroutine calculates the 
space-charge acceleration of the disk by summing the accelerations due to all 
the other disks. The acceleration of one disk due to another Is obtained by 
linear Interpolation on the space-charge force tables. The coefficients of 
the new quadratic curve fits are then calculated. Finally, the subroutine 
determines whether the frequency of the space-charge calculation should be 
changed. This Is done by comparing the space-charge accelerations based on 
the old curve fits with those just calculated. If the difference for any disk 
Is greater than the product of TOLSC and the maximum acceleration due to the 
rf electric field In the cavity, the frequency of calculating space-charge 
accelerations Is doubled. 


Subroutine SCFT3 

Subroutine SCFT3 Is called by subroutine ACC and In some cases by sub- 
routines INIT and TRN13. Subroutine SCFT3 Is called whenever a new calcula- 
tion of space-charge accelerations In the three-dimensional mode Is required. 
This subroutine first calculates what the space-charge accelerations would be 
If the present curve fits are used. Then for each ring the subroutine calcu- 
lates the space-charge acceleration of the ring by summing the accelerations 
due to all the other rings. The acceleration of one ring due to another Is 
obtained by three-dimensional linear Interpolation on the space-charge force 
tables. The coefficients of the new quadratic curve fits are then calculated. 
Finally, the subroutine determines whether the frequency of the space-charge 
calculation should be changed. This Is done by comparing the space-charge 
accelerations based on the old curve fits with those Just calculated. If the 
difference for any ring Is greater than the product of TOLSC and the maximum 
acceleration due to the rf electric field In the cavity, the frequency of cal- 
culating space-charge accelerations Is doubled. 


Subroutine SHSIG 

Subroutine SMSIG Is called by the main program. Subroutine SMSIG calcu- 
lates small-signal parameters as given In reference 7. This subroutine Is 
Identical to the one In the one-dimensional program described In reference 1 . 
A detailed description of the subroutine Is given In that reference. 


Subroutine STOAT 

Subroutine STDAT Is called by the main program. Subroutine STOAT calcu- 
lates standard data, that Is, constants whose values are Independent of Input 
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Subroutine STOB 


Subroutine STOB Is called by subroutine INDAT when the Initial cavity 
(INITCV) Is not one and when the Initial state Is not read from a dataset. 
Subroutine STOB obtains the Initial state from COMMON/STATE/. 


Subroutine TABLl 

Subroutine TABLl Is called In the one-dimensional mode by subroutine 
BTWNC when a new cavity is entered and tables have to be computed for the new 
cavity. If the new cavity has the same properties as the old cavity, sub- 
routine TABLl Is not called. This subroutine Is essentially the same as 
subroutine TABLE In the one-dimensional program described in reference 4. 

A detailed description of this subroutine Is given In that reference. 


Subroutine TBL3E 

Subroutine TBL3E Is called In the three-dimensional mode by subroutine 
BTWNC when a new cavity Is entered. Subroutine TBL3E computes tables for 
electric-field shapes and voltage-jump shapes. Also computed are the tables 
needed to calculate Induced voltages. If the new cavity has the same proper- 
ties as the old cavity, new tables are not required and subroutine TBL3E makes 
a return to subroutine BTWNC. 


Subroutine TBL3E calculates the following tables: 


TPXR(j,n,l) Real and Imaginary parts of the table Qz(n,z n ) for the shaping 
TPXI(j,n,l) factor of the axial electric field. See equation (62) for the 
definition of Q z (rj,z n ). 


TPXR(j,n,3) Real and Imaginary parts of the table Q r (rj,z n ) for the shaping 
TPXI(j,n,3) factor of the radial electric field. See equation (63) for the 
definition of Q r (rj,z n ). 


TPXR(j,n,2) Real and Imaginary parts of the table S(ri,z n ) used In calcu- 

TPXI(j,n,2) latlng Induced voltages. See equation (121) for the definition of 

s ( r j» z n)* 


TPXR(j,n,4) Table R z (rj,z n ) for the shaping factor of the axial voltage- 

jump field. See equation (76) for the definition of Rz^ r j» z n)- 

TPXR(j,n,5) Table Rr.(rj,z n ) for the shaping factor of the radial voltage- 
jump field. See equation (77) for the definition of R r (rj,z n ). 


Subroutine TBL3M 

Subroutine TBL3M Is used In the three-dimensional mode only. This sub- 
routine Is called by subroutine BTWNC when a new cavity Is entered and Is 
called by subroutine BTWNS when a new magnetic section Is entered. For single 
period PPM focusing, a magnetic section Is defined to be one-half of the 
magnetic period. For double period PPM focusing, a magnetic section Is de- 
fined to be one-fourth of the magnetic period. 
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Subroutine TBL3M calculates the following tables: 

BSUBZB(J,n) Axial magnetic field B z at the grid point (n,z n ) 

B$UBRB(j,n) Radial magnetic field B r at the grid point (n,z n ) 

PSIB(J,n) Magnetic flux field t at the grid point (rj,zj,) 

The axial grid points are equally spaced over the length of the magnetic 
section. If each magnetic section does not coincide with a cavity, these axial 
grid points will not coincide with the Integration nodes. Subroutine MAGCF Is 
called to obtain B z , B r , and »|» at the Integration nodes. 


Subroutine TRN13 

Subroutine TRN13 Is called by subroutine BTWNC when the transition from 
the one-dimensional mode to the three-dimensional mode Is made. Subroutine 
TRN13 recalulates all those variables whose values change as a result of this 
transition. 


Subroutine XINT 

Subroutine XINT Is called by the main program. Subroutine XINT Integrates 
the equations of motion over one step. The numerical Integration method Is 
described In the section on the dynamics of beam rings. 


CONCLUDING REMARKS 

The use of the coupled-cavity traveling wave tube for space communica- 
tions had led to an Increased Interest In Improving the efficiency of the 
basic Interaction process In these devices through velocity resynchronization 
and other methods. In order to analyze these methods, we have recently 
developed a flexible, three-dimensional, axlally-symmetrlc, large-signal 
computer program for use on the IBM 370 time-sharing system. The present 
report Is a users' manual for this program. The report describes the program 
In sufficient detail to allow a user to make modifications In the program If 
desired. 
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APPENDIX A 
SYMBOLS 


a tunnel radius 


j 


a 

a. 


b 

c 



(r,z) 


m 


N 

N 

N 

P 


c 

d 

z 


Q z (r,z) 

Q r (r.z) 

q 

R 

R z (r,z) 

R r (r,z) 

r(z) 

r 1n 

r out 

Ar 

sc 

S k (r,z) 

t(z) 

T 


Inner diameter of ring magnet 

n *^ 1 Fourier coefficient In expansion for magnetic field 

magnitude of axial magnetic field at r = a and z = middle of gap 

radial component of magnetic field 

axial component of magnetic field 

Initial beam radius 

velocity of light 

half-length of magnetic gap 

modified Bessel function of order s 

fundamental Fourier component of the beam current density 

Bessel function of order s 

length of cavity 

axial thickness of a ring 

half-length of k ** 1 gap 

mass of a ring 

number of cavities 

number of disks per beam wavelength 

number of nodes per cavity 

magnetic period 

table for calculation of rf electric field forces 

table for calculation of rf electric field forces 

charge of a ring 

number of rings per disk 

table for calculation of voltage-jump forces 

table for calculation of voltage-jump forces 

centroid radius of a ring at axial location z 

Inner radius of a ring 

outer radius of a ring 

radial grid spacing for space charge tables 
table for calculation of Induced voltages 
time of arrival of a ring at axial location z 
reciprocal of frequency 
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u q Initial beam velocity 

V Q beam voltage 

V k gap voltage, k th cavity 

V bk backward voltage, k th cavity 

V fk forward voltage, k th cavity 

AV bk backward Induced voltage, k th cavity 

AV fk forward Induced voltage, k th cavity 

Zk Interaction Impedance for k th cavity 

2 axial velocity of a ring 

^ z sc axial grid spacing for space charge tables 

Az k Integration step size In k th cavity 

z nk ax1al Position, relative to beginning of first cavity, of n th node 

In k** 1 cavity 

(<*L) b k loss factor for backward voltage In k th cavity 

(°>L)fk loss factor for forward voltage In k th cavity 

propagation factor, defined In equation ( 60 ) 

( B-j L) k phase shift of voltage for k** 1 cavity 
permittivity of free space 

eU) normalized time of arrival of a ring at normalized axial location 5 

beam wavelength 

v shaping parameter for electric field 

v shaping factor for magnetic field 

l normalized axial position; Independent variable for equations of 

motion 

pU) normalized centroid radius of a ring at normalized axial location 5 

<j(z) charge density of a ring at axial location z 

<t> angular velocity of a ring 

V magnetic flux 

u angular frequency 
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APPENDIX B 


DERIVATION OF FOURIER COMPONENT OF BEAM CURRENT DENSITY 

In a previous section, the fundamental Fourier component of the beam cur- 
rent density was defined as 


J l7 (r,z) = i J (r,z,f)e' ,ut ' dt * (B-l) 

,z 1 t-T z 

We proceed to obtain an expression for J lz (r,z) which Is convenient for 
computation. The following analysis Is a simple generalization of the one- 
dimensional analysis given In appendix B of reference 5. Let tT(z) be the 

1L * 

time of arrival at z of the 1 ring In the nr radiofrequency cycle. 

We let m = 0 for the radiofrequency cycle considered In the program. To con- 
form with previous notation, we will henceforth omit the 0 subscript when 

referring to the cycle considered In the program; that Is, t°(z) = t^(z). 

Then for any Integer m. 


t"(z) = t,(z) - mT 


( B-2) 


In the following discussion we will use the rectangular function 
g(x; X] ,X 2 ) defined as 


g(x; x r x 2 ) = 1, 


x ] < x < x^ 


g(x; x r x 2 ) = 0, 


x < x 1 or x > x ? 


(B-3) 


Now consider J^(r,z,t) as a function of t for some fixed r and z 
(see fig. 12). We can visualize this as an observer at (r,z), observing 
rings as they pass by. At t = t^z), the value of J z (r,z,t) Is 
«,<z)v,z(z)g< r ; r^ ^.r^ out ). where o^z) Is the charge density of the 1 th 
ring. v^(z) the axial velocity of the I** 1 ring, and r^ and r^ 

are the Inner and outer radius of the I** 1 ring. Furthermore, J z (r,z,t) has 
this value for a time Interval centered about t^(z), where the length of the 
Interval Is l^/v^fz). We can express this partial contribution to 
J z (r,z,t) as 


o 1 (z)v 1z (z)g(r; r^.r^g^t; t*(z) - t "(z) ♦ 27^(1)) < B " 4 > 

The total contribution to J z (r,z,t) Is obtained by summing over all rings 
and cycles. 
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N d R 


^ W u lzMn r 1j1n ,r 10ut ) 


1=1 m = -■ 


x g 


( t; t: 


m (Z ) _ — t m rzi + * d \ 

'' ’ 2v 1z (z)« V z) + 2v 1z (z)^ 


We the have for Jiz(r,z) 

N d R 

V r - Z > - T L Z °1< z >V z >S< r : r l,ln> r i.out ) 


1=1 m = 


(B-5) 


/ g^f; t™ 


<z) -S^)' *?**> 


) 


( B-6) 


By the perlodlcty of the motion, there will be N d R rings passing the posi- 
tion z In the time Interval from t-T to t. In general, the N d R rings 
do not all come from the same rf cycle. Let *m^ be the cycle number of the 
1 tn ring. From equations (B-2) and (B-6), we obtain 


N d R 

’u (r ' Z) ■ T E o 1 (Ov lz (z)g(r; r, , n .r , t ) 

1=1 *11 


12 .- Wt ' df (B-7) 


where t^j and t ^2 are given by 


*11 = *1 (Z) “ m 1 T " 2v 1z (z) 


( B— 8) 


t^ = t^(z) - m,T + 


1 2v 1z (z) 


Evaluating the Integral In equation (B-7) yields 


N d R 


,J lz (r,2) ■ T 2 °1 (2)ll d 9(r; out) 


sin 


l \ 

\ 2v 1z< 2 )/ 


( B— 9 ) 


2v 1z (z)y -Mt^j-mj] 

“ l d X 6 


1 = 1 


2v 1z< z > 


( B— 1 0) 
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Since com^T = 2irm^, we have 


j, z < r 


, > \2v^ {z)J -lut,(z) 

• z > ■ T < ’l (z> ‘d s<ri r 1,1n- r 1,out> ST * 6 


1 = 1 


?v u (z) 


The charge density o\(z) is given by 


o,(z) = / 2 2 \ 

*( r 1,out ” r 1,1nj l d 

where q^(z) Is the charge of the 1^ ring. Combining equations (B-11) 
and 

(B-12) yields 


q^z) 


J ' 2<r, ° = f S Jr 2 . 


9(r; r 1.1n- r 1.out )1 'i (z) 

(2 _2 \ 


1 = 1 


2 2 \ 
1,out " r 1,1nj 


\2v 1z (z)/ -Iwt^z) 
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Figure 3. - Nodes and antinodes for cavity k. 



Axial space-charge force on one Innermost 
due to another 



Figure 4 - Axial space-charge force on one innermost ring due to another 
as a function of normalized axial separation distance; tJ a ■ 1/3. 


Comparison of R-tables model to R ^-tables model of space charge forces. 
Consider 3 rings, all at same axial position, and all having their initial 
dimensions. Beam radius ■ a 7 x tube radius. Innermost ring is ring 
#1, etc. Compare radial space charge forces. Axial thickness is y 24 
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Figure 5. * Typical comparison of R 2 -tables model to R-tables model. 
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Figure 9. - Trajectory obtained from 
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Figure 10. - Comparison of trajectory obtained from exact solution of equation r" ■ K/r to the trajectory 
obtained from model used in computer program. 
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